Statistical Mechanics and the Physics of the 
Many-Particle Model Systems* 



The development of methods of quantum statistical mechanics is considered in light of 
their applications to quantum solid-state theory. We discuss fundamental problems of the 
physics of magnetic materials and the methods of the quantum theory of magnetism, includ- 
ing the method of two-time temperature Green's functions, which is widely used in various 
physical problems of many-particle systems with interaction. Quantum cooperative effects and 
quasiparticle dynamics in the basic microscopic models of quantum theory of magnetism: the 
Heisenberg model, the Hubbard model, the Anderson Model, and the spin-fermion model are 
considered in the framework of novel self-consistent-field approximation. We present a com- 
parative analysis of these models; in particular, we compare their applicability for description 
of complex magnetic materials. The concepts of broken symmetry, quantum protectorate, 
and quasiaverages are analyzed in the context of quantum theory of magnetism and theory 
of superconductivity. The notion of broken symmetry is presented within the nonequilibrium 
statistical operator approach developed by D.N. Zubarev. In the framework of the latter ap- 
proach we discuss the derivation of kinetic equations for a system in a thermal bath. Finally, 
the results of investigation of the dynamic behavior of a particle in an environment, taking 
into account dissipative effects, are presented. 
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1 Introduction 



The purpose of this review is to trace the development of some methods of quantum statistical 
mechanics formulated by N.N. Bogoliubov, and also to show their effectiveness in applications to 
problems of quantum solid-state theory, and especially to problems of quantum theory of mag- 
netism. It is necessary to stress, that the path to understanding the foundations of the modern 
statistical mechanics and the development of efficient methods for computing different physical 
characteristics of many-particle systems was quite complex. The main postulates of the mod- 
ern statistical mechanics were formulated in the papers by J. P. Joule (1818-1889), R. Clausius 
(1822-1888), W. Thomson (1824-1907), J.C. Maxwell (1831-1879), L. Boltzmann (1844-1906), 
and, especially, by J.W. Gibbs (1839-1903). The monograph by Gibbs "Elementary Principles 
in Statistical Mechanics Developed with Special Reference to the Rational Foundations of Ther- 
modynamics"^^ remains one of the highest peaks of modern theoretical science. A significant 
contribution to the development of modern methods of equilibrium and nonequilibrium statistical 
mechanics was made by Academician N.N. Bogoliubov (1909-1992) 

Specialists in theoretical physics, as well as experimentalists, must be able to find their way 
through theoretical problems of the modern physics of many-particle systems because of the fol- 
lowing reasons. Firstly, the statistical mechanics is filled with concepts, which widen the physical 
horizon and general world outlook. Secondly, statistical mechanics and, especially, quantum statis- 
tical mechanics demonstrate remarkable efficiency and predictive ability achieved by constructing 
and applying fairly simple (and at times even crude) many-particle models. Quite surprisingly, 
these simplified models allow one to describe a wide diversity of real substances, materials, and 
the most nontrivial many-particle systems, such as quark-gluon plasma, the DNA molecule, and 
interstellar matter. In systems of many interacting particles an important role is played by the 
so-called correlation effects,^ which determine specific features in the behavior of most diverse 
objects, from cosmic systems to atomic nuclei. This is especially true in the case of solid-state 
physics. Investigation of systems with strong inter-electron correlations, complicated character of 
quasiparticle states, and strong potential scattering is an extremely important and topical problem 
of the modern theory of condensed matter. Our time is marked by a rapid advancement in design 
and application of new materials, which not only find a wide range of applications in different 
areas of engineering, but they are also connected with the most fundamental problems in physics, 
physical chemistry, molecular biology, and other branches of science. The quantum cooperative 
effects, such as magnetism and superconductivity, frequently determine the unusual properties of 
these new materials. The same can be also said about other non-trivial quantum effects like, for 
instance, the quantum Hall effect, the Bose-Einstein condensation, quantum tunneling and oth- 
ers. This research direction is developing very rapidly, setting a fast pace for widening the domain 
where the methods of quantum statistical mechanics are applied. This review will support the 
above statement by concrete examples. 

2 Quantum Statistical Mechanics and Solid State Physics 

The development of experimental techniques over the recent years opened the possibility for 
synthesis and investigations of a wide class of new substances with unusual combination of prop- 
erty EHE] Transition and rare-earth metals and especially compounds containing transition and 
rare-earth elements possess a fairly diverse range of properties. Among those, one can men- 
tion magnetically ordered crystals, superconductors, compounds with variable valence and heavy 
fermions, as well as substances which under certain conditions undergo a metal-insulator transi- 
tion, like perovskite-type manganites, which possesses a large magneto-resistance with a negative 
sign. These properties find widest applications in engineering; therefore, investigations of this 
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class of substances should be classified as the currently most important problems in the physics 
of condensed matter. A comprehensive description of materials and their properties (as well as 
efficient predictions of properties of new materials) is possible only in those cases, when there is an 
adequate quantum-statistical theory based on the information about the electron and crystalline 
structures. The main theoretical problem of this direction of research, which is the essence of 
the quantum theory of magnetismj^^ is investigations and improvements of quantum-statistical 
models describing the behavior of the above-mentioned compounds in order to take into account 
the main features of their electronic structure, namely, their dual "band-atomic" natureP^E^ 
The construction of a consistent theory explaining the electronic structure of these substances 
encounters serious difficulties when trying to describe the collectivization- localization duality in 
the behavior of electrons. This problem appears to be extremely important, since its solution 
gives us a key to understanding magnetic, electronic, and other properties of this diverse group of 
substances. The author of the present review investigated the suitability of the basic models with 
strong electron correlations and with a complex spectrum for an adequate and correct description 
of the dual character of electron states. A universal mathematical formalism was developed for 
this investigation.^ It takes into account the main features of the electronic structure and allows 
one to describe the true quasiparticle spectrum, as well as the appearance of the magnetically or- 
dered, superconducting, and dielectric (or semiconducting) states. With a few exceptions, diverse 
physical phenomena observed in compounds and alloys of transition and rare- earth metals ™EB 
cannot be explained in the framework of the mean-field approximation, which overestimates the 
role of inter-electron correlations in computations of their static and dynamic characteristics. The 
circle of questions without a precise and definitive answer, so far, includes such extremely impor- 
tant (not only from a theoretical, but also from a practical point of view) problems as the adequate 
description of quasiparticle dynamics for quantum-statistical models in a wide range of their pa- 
rameter values. The source of difficulties here lies not only in the complexity of calculations of 
certain dynamic properties (such as, the density of states, electrical conductivity, susceptibility, 
electron-phonon spectral function, the inelastic scattering cross section for slow neutrons), but also 
in the absence of a well-developed method for a consistent quantum-statistical analysis of a many- 
particle interaction in such systems. A self-consistent field approach was used in the papers^^HEZl 
for description of various dynamic characteristics of strongly correlated electronic systems. It al- 
lows one to consistently and quite compactly compute quasiparticle spectra for many-particle 
systems with strong interaction taking into account damping effects. The correlation effects and 
quasiparticle damping are the determining factors in analysis of the normal properties of high- 
temperature superconductors, and of the transition mechanism into the superconducting phase. 
We also formulated a general scheme for a theoretical description of electronic properties of many- 
particle systems taking into account strong inter-electron correlations] 20 ^ 22 " 24 ^ 26 ^ 2 ' The scheme is 
a synthesis of the method of two-time temperature Green's functions^ and the diagram technique. 
An important feature of this approach is a clear-cut separation of the elastic and inelastic scatter- 
ing processes in many-particle systems (which is a highly nontrivial task for strongly correlated 
systems). As a result, one can construct a correct basic approximation in terms of generalized 
mean fields (the elastic scattering corrections), which allows one to describe magnetically ordered 
or superconducting states of the system. The residual correlation effects, which are the source of 
quasiparticle damping, are described in terms of the Dyson equation with a formally exact repre- 
sentation for the mass operator. There is a general agreement that for heavy-fermion compounds, 
the model Hamiltonian is well established (the periodic Anderson model or the periodic Kondo 
lattice), and the main theoretical challenge in this case lies in constructing accurate approxima- 
tions. However, in the case of high-temperature superconductors or perovskite-type manganites, 
neither a model, nor adequate approximate analytical methods for its solution are available. Thus, 
the development and improvement of the methods of quantum statistical mechanics still remains 
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quite an important direction of researchP^l 

3 Magnetic Properties of Substances and Models of Magnetic 
Materials 

It is widely accepted that the appearance of magnetically ordered states in transition metals is 
to some extent a consequence of the atom-like character of d-states, but mostly it is the result of 
interatomic exchange interactions. In order to better understand the origin of quantum models of 
magnetic materials, we discuss here briefly the physical aspects of the magnetic properties of solid 
materials. The magnetic properties of substances belong to the class of natural phenomena, which 
were noticed a long time agoP^^ Although it is assumed that we encounter magnetic natural 
phenomena less frequently than the electric ones, nevertheless as was noticed by V. Weisskopf, 
"magnetism is a striking phenomenon; when we hold a magnet in one hand and a piece of iron 
in another, we feel a peculiar force, some "force of nature", similar to the force of gravity".^ 
It is interesting to note that the experiment-based scientific approach began from investigations 
of magnetic materials. This is the so called inductive method, which insists on searching for 
truth about the nature not in deductions, not in syllogisms and formal logics, but in experiments 
with the natural substances themselves. This method was applied for the first time by William 
Gilbert (1544-1603), Queen Elizabeth's physician. In his book, "On the magnet, magnetic bodies, 
and on the great magnet, the Earth" published in 1600, he described over 600 specially per- 
formed experiments with magnetic materials, which had led him to an extremely important and 
unexpected for the time conclusion, that the Earth is a giant spherical magnet. Investigations 
of Earth's and other planet's magnetism is still an interesting and quite important problem of 
modern scienceP^HHH Thus, it is with investigations of the physics of magnetic phenomena, that 
the modern experiment based science began. Note, that although the creation of the modern 
scientific methods is often attributed to Francis Bacon, Gilbert's book appeared 20 years earlier 
than "The New Organon" by Francis Bacon (1561-1626). 

The key to understanding the nature of magnetism became the discovery of a close connection be- 
tween magnetism and electricity.^ For a long time the understanding of magnetism's nature was 
based on the hypothesis of how the magnetic force is created by magnets. Andre Ampere (1775- 
1836) conjectured that the principle behind the operation of an ordinary steel magnet should be 
similar to an electric current passing over a circular or spiral wire. The essence of his hypothesis 
laid in the assumption, that each atom contains a weak circular current, and if most of these 
atomic currents are oriented in the same direction, then the magnetic force appears. All subse- 
quent developments of the theory of magnetism consisted in the development and refinement of 
Ampere's molecular currents hypothesis. 

As an extension of this idea by Ampere a conjecture was put forward that a magnet is an ensem- 
ble of elementary double poles, magnetic dipoles. The dipoles have two magnetic poles which are 
inseparably linked. In 1907, Pierre Weiss (1865-1940) proposed a phenomenological picture of the 
magnetically ordered state of matter. He was the first to perform a phenomenological quantita- 
tive analysis of the magnetic phenomena in substances.^ Weiss's investigations were based on 
the notion, introduced by him, of a molecular field. Subsequently this approach was named the 
molecular (or mean, or effective) field approximation, and it is still being widely used even at the 
present time.^ The simplest microscopic model of a ferromagnet in the molecular field's approx- 
imation is based on the assumption that electrons form a free gas of magnetic arrows (magnetic 
dipoles), which imitate Ampere's molecular currents. In the simplest case it is assumed that these 
"elementary magnets" could orient in space either along a particular direction, or against it. In 
order to find the thermodynamically equilibrium value of the magnetization (M) as a function of 
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temperature T, one has to turn to general laws of thermodynamics. This is especially important 
when we consider the behavior of a system at finite temperatures. Finding the equilibrium mag- 
netization of a ferromagnet becomes quite a simple task if we first succeed in writing down its 
energy E({M)) as a function of magnetization. All we have to do after that is to minimize the 
free energy F((M)), which is defined by the following relationship!^ 

F{(M)) = E((M)) -TS((M)). (3.1) 

Here, S({M)) is the system entropy also written down as a function of magnetization. It is im- 
portant to stress, that the problem of calculating the system entropy cannot be solved in the 
framework of just thermodynamics. In order to find the entropy one has to turn to statistical 
mechanics jHIMEl which provides a microscopic foundation to the laws of thermodynamics. Note 
that derivations of equilibrium magnetization (M) as a function of temperature T, or, more gener- 
ally, investigations of relationships between the free energy and the order parameter in magnetics 
and pyroelectrics are still ongoing even at the present time! 42 " 45 ! Of course, modern investigations 
take into account all the previously accumulated experience. In the framework of the P. Weiss 
approach one investigates the appearance of a spontaneous magnetization (M) ^ for H = 0. 
This approach is based on the following postulate for the behavior of the system's energy 

E({M)) ~NI({M)) 2 . (3.2) 

This expression takes into account the interaction between elementary magnets (arrows). Here, 
I is the energy of the Weiss molecular field per atomic magnetic arrow. The question on the 
microscopic nature of this field is beyond the framework of the Weiss approach. The minimization 
of the free energy F((M)) yields the following relationship 

(M) = th(T c (M)/T), (3.3) 

Where Tq is the Curie or the critical temperature. As the temperature decreases below this 
critical value a spontaneous magnetization appears in the system. The Curie temperature was 
named in honor of Pierre Curie (1859-1906), who established the following law for the behavior 
of susceptibility x i n paramagnetic substances 

X = lim <^> = %. (3.4) 

Depending on the actual material, the Curie parameter C obtains different (positive) values.^ 
Note that Pierre Curie performed thorough investigations of the magnetic properties of iron back 
in 1895. In the process of those experiments he established the existence of a critical temper- 
ature for iron, above which the ferromagnetic properties disappear. These investigations laid a 
foundation for investigations of order-disorder phase transitions, and other phase transformations 
in gases, liquids, and solid substances. This research direction created the core of the physics of 
critical phenomena, which studies the behavior of substances in the vicinity of critical tempera- 
tures. S3 

Extensive investigations of spontaneous magnetization and other thermal effects in nickel in the 
vicinity of the Curie temperature were performed by Weiss and his collaborators.^ They devel- 
oped a technique for measuring the behavior of the spontaneous magnetization in experimental 
samples for different values of temperature. Knowing the behavior of spontaneous magnetization 
as a function of temperature one can determine the character of magnetic transformations in the 
material under investigation. Investigations of the behavior of the magnetic susceptibility as a 
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function of temperature in various substances remain important even at the present timepSHSS 
Within the P. Weiss approach we obtain the following expression for the Curie temperature 

T c = 21 /k B . (3.5) 

In order to obtain a rough estimate for the magnitude of / take Tc = 1000 K. Then one obtains 
/ ~ 10 -13 erg/atom. This implies that the only origin of the Weiss molecular field can be the 
Coulomb interaction of electrical charges P^"^ Computations in the framework of the molecular 
field method lead to the following formula for the magnetic susceptibility 

X - N A (M)/H _ 1 -^ m , (3.6) 

where [is = eh/2mc is the Bohr magneton (within the Weiss approach this is the magnetic 
moment of the magnetic arrows imitating Ampere's molecular currents). The above expression 
for the susceptibility is referred to as the Curie- Weiss law. Thus, the Weiss molecular field, whose 
magnitude is proportional to the magnetization, is given by 

H w = k B T c (M)/n B . (3.7) 

Researchers tried to find the answer to the question on the nature of this internal molecular field 
in ferromagnets for a long time. That is, they tried to figure out which interaction causes the 
parallel alignment of electron spins. As was stressed in the book: 51 "At first researchers tried to 
imagine this interaction of electrons in a given atom with surrounding electrons as some quasi- 
magnetic molecular field, acting on the electrons of the given atom. This hypothesis served as a 
foundation for the P. Weiss theory, which allowed one to describe qualitatively the main properties 
of ferromagnets". However, it was established that Weiss's molecular field approximation is appli- 
cable neither for theoretical interpretation, nor for quantitative description of various phenomena 
taking place in the vicinity of the Curie temperature. Although numerous attempts aiming to 
improve Weiss's mean-field theory were undertaken, none of them led to significant progress in 
this direction. Numerical estimates yield the value H\y = 10 7 oersted for the magnitude of the 
Weiss mean field. The nonmagnetic nature of the Weiss molecular field was established by direct 
experiments in 1927 (see the books^E}. Ya.G. Dorfman performed the following experiment. 
An electron beam passing though nickel foil magnetized to the saturation level is falling on a 
photographic plate. It was expected that if such a strong magnetic field indeed exists in nickel, 
then the electrons passing through the magnetized foil would deflect. However, it turned out that 
the observed electron deflection is extremely small. The experiment led to the conclusion that, 
contrary to the consequences of the Weiss theory, the internal fields of large intensity are not 
present in ferromagnets. Therefore, the spin ordering in ferromagnets is caused by forces of a 
nonmagnetic origin. It is interesting that fairly recently, in 2001, similar experiments were per- 
formed again 52 (in a substantially modified form, of course). A beam of polarized "hot" electrons 
was scattered by thin ferromagnetic nickel, iron, and cobalt films, and the polarization of the 
scattered electrons was measured. The concept of the Weiss exchange field W(x) ~ — J a S(x) was 
used for theoretical analysis.^E^ The real part of this field corresponds to the exchange inter- 
action between the incoming electrons and the electron density of the film (the imaginary part 
is responsible for absorption processes). The derived equations, describing the beam scattering, 
resemble quite closely the corresponding equations for the Faraday's rotation effect in the light 
passing through a magnetized environment ! 52 H 53 I The theoretical consideration is based on using 
the mean-field approximation, namely on the replacement 

W ~ (W(a?)) = J a (S{x)). (3.8) 
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The subsequent quite rigorous and detailed consideration^ aimed at deriving the effective quan- 
tum dynamics of the field W(x) showed, that this dynamics is described by the Landau-Lifshitz 
equation.^ The spatial and temporal variations of the field W(x) are described by spin waves. 
The quanta of the Weiss exchange field are magnons. 

One has to note that in its original version, the Weiss molecular field was assumed to be uniformly 
distributed over the entire volume of the sample, and had the same magnitude in all points of 
the substance. An entirely different situation takes place in a special class of substances called 
antiferromagnets. As the temperature of antiferromagnets falls below a particular value, a mag- 
netically ordered state appears in the form of two inserted into each other sublattices with opposite 
directions of the magnetization. This special value of the temperature became known as the Neel 
temperature, after the founder of the antiferromagnetism theory L. Neel (1904-2000). In order to 
explain the nature of the antiferromagnetism (as well as of the ferrimagnetism) L. Neel introduced 
a profound and nontrivial notion of local molecular fields. 54 However, there was no a unified ap- 
proach to investigations of magnetic transformations in real substances. Moreover, a consistent 
consideration of various aspects of the physics of magnetic phenomena on the basis of quantum 
mechanics and statistical physics was and still is an exceptionally difficult task, which to the 
present days does not have a complete solutionP^E^ This was the reason why the authors of the 
most complete, at that time, monograph on the magnetism characterized the state of affairs in the 
physics of magnetic phenomena as follows: "Even recently the problems of magnetism seemed to 
belong to an exceptionally unrewarding area for theoretical investigations. Such a situation could 
be explained by the fact that the attention of researchers was devoted mostly to ferromagnetic 
phenomena, because they played and still play quite an important role in engineering. However, 
the theoretical interpretation of the ferromagnetism presents such formidable difficulties, that to 
the present day this area remains one of the darkest spots in the entire domain of physics".^ 
The magnetic properties and the structure of matter turned out to be interconnected subjects. 
Therefore, a systematic quantum-mechanical examination of the problem of magnetic substances 
was considered by most reseaxcherPB^HSQ] as q U ite an important task. Heisenberg, Dirac, Hund, 
Pauli, van Vleck, Slater, and many other researchers contributed to the development of the quan- 
tum theory of magnetism. As was noted by D. MattisJ 13 "...by 1930, after four years of most 
exciting and striking discoveries in the history of theoretical physics, a foundation for the modern 
electron theory of matter was laid down, after that an epoch of consolidation and computations 
had began, which continues up to the present day" . 

Over the last decades the physics of magnetic phenomena became a vast and ramified domain of 
modern physical science.^ ^^ss M^ 4 ^ The rapid development of the physics of magnetic materi- 
als was influenced by introduction and development of new physical methods for investigating the 
structural and dynamical properties of magnetic substances.^ These methods include magnetic 
neutron diffraction analysisj^"^ NMR and EPR-spectroscopy, the Mossbauer effect, novel opti- 
cal methods,^ as well as recent applications of synchrotron radiation.^HEZl T n particular, unique 
possibilities of the thermal neutron's scattering methods^^^Sl allow one to obtain information 
on the magnetic and crystalline structure of substances, on the distribution of magnetic moments, 
on the spectrum of magnetic excitations, on critical fluctuations, and on many other properties 
of magnetic materials. In order to interpret the data obtained via inelastic scattering of slow 
neutrons one has to take into account electron-electron and electron-nuclear interactions in the 
system, as well as the Pauli exclusion principle. Here, we again face the challenge of consider- 
ing various aspects of the physics of magnetic phenomena, consistently on the basis of quantum 
mechanics and statistical physics. In other words, we are dealing with constructing a consistent 
quantum theory of magnetic substances. As was rightly noticed by K. Yosida, "The question of 
electron correlations in complex electronic systems is the beginning and the end of all research on 
magnetism" Thus, the phenomena of magnetism can be described and interpreted consistently 
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only in the framework of quantum statistical theory of many interacting particles. 

4 Quantum Theory of Magnetism 

It is well known that "quantum mechanics is the key to understanding magnetism" One of 
the first steps in this direction was the formulation of "Hund's rules" in atomic physicsP^ As 
was noticed by D. Mattis, 17 "The accumulated spectroscopic data allowed Stoner (1899-1968) to 
attribute the correct number of equivalent electrons to each atomic shell, and Hund (1896-1997) 
to state his rules, related to the spontaneous magnetic moments of a free atom or ion". Hund's 
rules are empirical recipes. Their consistent derivation is a difficult task. These rules are stated 
as follows: 

(1) The ground state of an atom or an ion with a L — S coupling is a state with the maximal 
multiplicity (2S + 1) for a given electron configuration. 

(2) From all possible states with the maximal multiplicity, the ground state is a state with the 
maximal value of L allowed by the Pauli exclusion principle. 

Note that the applicability of these empirical rules is not restricted to the case, when all electrons 
lie in a single unfilled valency shell. A rigorous derivation of Hund's rules is still missing. How- 
ever, there are a few particular cases which show their validity under certain restrictions^ 3 -"^ (see 
a recent detailed analysis of this question in the paperJ^"^). Nevertheless Hund's rules are very 
useful and are widely used for analysis of various magnetic phenomena. A physical analysis of 
the first Hund's rule leads us to the conclusion, that it is based on the fact, that the elements of 
the diagonal matrix of the electron-electron's Coulomb interaction contain the exchange's inter- 
action terms, which are entirely negative. This is the case only for electrons with parallel spins. 
Therefore, the more electrons with parallel spins involved, the greater the negative contribution 
of the exchange to the diagonal elements of the energy matrix. Thus, the first Hund's rule implies 
that electrons with parallel spins "tend to avoid each other" spatially. Here, we have a direct 
connection between Hund's rules and the Pauli exclusion principle. 

One can say that the Pauli exclusion principle (1925) lies in the foundation of the quantum theory 
of magnetic phenomena. Although this principle is merely an empirical rule, it has deep and 
important implications.^ W. Pauli (1900-1958) was puzzled by the results of the ortho-helium 
terms analysis, namely, by the absence in the term structure of the presumed ground state, that 
is the (l 3 ^) level. This observation stimulated him to perform a general examination of atomic 
spectra, with the aim to find out, if certain terms are absent in other chemical elements and under 
other conditions as well. It turned out, that this was indeed the case. Moreover, the conducted 
analysis of term systems had shown that in all the instances of missing terms the entire sets of 
the quantum numbers were identical for some electrons. And vice versa, it turned out that terms 
always drop out in the cases when entire sets of quantum numbers are identical. This observation 
became the essence of the Pauli exclusion principle: 

The sets of quantum numbers for two (or many) electrons are never identical; two sets of quantum 
numbers, which can be obtained from one another by permutations of two electrons, define the 
same state. 

In the language of many-electron wave functions one has to consider permutations of spatial and 
spin coordinates of electrons i and j in the case when both the spin variables o~i = o~j = o~q and 
the spatial coordinates r*j = rj = r*o of these two electrons are identical. Then, we obtain: 

Pijip(fiai, . . . no-i, . . . rjo-j, ■■■) = ip{nai, . . . noi, . . . rjo-j, . . .). (4.1) 

The Pauli exclusion principle implies that 

Pijipinai, . . . fi<T i} . . . rjo-j, ■■■) = -^(r>i, . . . r^, . . . fjaj, . . .). (4.2) 
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The above conditions are satisfied simultaneously only in the case, when ip is equal to zero iden- 
tically. Therefore, we arrive at the following conclusion: electrons are indistinguishable, that 
is, their permutations must not change observable properties of the system. The wave function 
changes or retains its sign under permutations of two particles depending on whether these indis- 
tinguishable particles are bosons or fermions. A consequence of the Pauli exclusion principle is 
the Aufbau principle, which leads to the periodicity in the properties of chemical elements. The 
fact that not more than one electron can occupy any single state leads also to such fundamental 
consequences as the very existence of solid bodies in nature. If the Pauli exclusion principle was 
not satisfied, no substance could ever be in a solid state. If the electrons would not have spin 
(that is, if they were bosons) all substances would occupy much smaller volumes (they would have 
higher densities), but they would not be rigid enough to have the properties of solid bodies. 
Thus, the tendency of electrons with parallel spins "to avoid each other" reduces the energy of 
electron-electron Coulomb interaction, and hence, lowers the system energy. This property has 
many important implications, in particular, the existence of magnetic substances. Due to the 
presence of an internal unfilled nd- or n/-shell, all free atoms of transition elements are strong 
magnetic, and this is a direct consequence of Hund's rules. When crystals are formed^E^'^'^ 
the electronic shells in atoms reorganize, and in order to understand clearly the properties of crys- 
talline substances, one has to know the wave function and the energies of (previously) outer-shell 
electrons. At the present time there are well-developed efficient methods for computing electronic 
energy levels in crystalsP^^ Speaking qualitatively, we have to find out how the atomic wave's 
functions change when crystals are formed, and how significantly they delocalizeP^ 



4.1 The Method of Model Hamiltonians 

The method of model Hamiltonians proved to be very efficient in the theory of magnetism. With- 
out any exaggeration one can say, that the tremendous successes in the physics of magnetic phe- 
nomena were achieved, largely, as a result of exploiting a few simple and schematic model concepts 
for " the theoretical interpretation of ferromagnetism" One can regard the Ising modeP^^l as 
the first model of the quantum theory of magnetism. In this model, formulated by W. Lenz (1888- 
1957) in 1920 and studied by E. Ising (1900-1998), it was assumed that the spins are arranged at 
the sites of a regular one-dimensional lattice. Each spin can obtain the values zLH/2: 

'H = -Y J hjSfS*. (4.3) 

<ij> 

This was one of the first attempts to describe the magnetism as a cooperative effect. It is interest- 
ing that the one-dimensional Ising model was solved by Ising in 1925, while the exact solution of 
the Ising model on a two-dimensional square lattice was obtained by L. Onsager ( 1903-1976 )(26J122I 
only in 1944. However, the Ising model oversimplifies the situation in real crystals. W. Heisen- 
berg (1901-1976)P^ and P. Dirac (1902-1984^2] formulated the Heisenberg model, describing the 
interaction between spins at different sites of a lattice by the following isotropic scalar function 

U = - Ai - j)SiSj - g^ B H Sf. (4.4) 

ij i 

Here J(i — j) (the "exchange integral") is the strength of the exchange interaction between the 
spins located at the lattice sites % and j. It is usually assumed that J(i — j) = J(j — i) and 
J(i ~ j = 0) = 0) which means that only the inter-site interaction is present (there is no self- 
interaction). The Heisenberg Hamiltonian (|4.4p can be rewritten in the following form: 

H = -Y t J{i- j)(SfS? + S+S7). (4.5) 

ij 
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Here, 5 = S x db iS y are the spin raising and lowering operators. They satisfy the following set 
of commutation relationships: 

[S+,Sr]_ = 2Sf6 ij ; [5+ 5~] + = 25(5 + 1) - 2(5f) 2 ; 
[S?, 5|]_ = ±Sf8 ir , Sf = 5(5 + 1) - (Sf) 2 - 5r5+; 

(5+) 2S+1 = 0, (5r) 25+1 = 0. 

Note that in the isotropic Heisenberg model the z-component of the total spin 5f ot = ^ Sf is a 
constant of motion, that is [H, Sf ot ] = 0. 

Thus, in the framework of the Heisenberg-Dirac-van Vleck modelpSEBHlou describing the interac- 
tion of localized spins, the necessary conditions for the existence of ferromagnetism involve the 
following two factors. Atoms of a "ferromagnet to be" must have a magnetic moment, arising 
due to unfilled electron d- or /-shells. The exchange integral Jjj related to the electron exchange 
between neighboring atoms must be positive. Upon fulfillment of these conditions the most en- 
ergetically favorable configurations in the absence of an external magnetic field correspond to 
parallel alignment of magnetic moments of atoms in small areas of the sample (domains) Of 
course, this simplified picture is only schematic. A detail derivation of the Heisenberg-Dirac-van 
Vleck model describing the interaction of localized spins is quite complicated. Because of a short- 
age of space we cannot enter into discussion of this quite interesting topic.EStUM' An important 
point to keep in mind here is that magnetic properties of substances are born by quantum effects, 
the forces of exchange interaction.^^ 

As was already mentioned above, the states with antiparallel alignment of neighboring atomic 
magnetic moments are realized in a fairly wide class of substances. As a rule, these are various 
compounds of transition and rare-earth elements, where the exchange integral Jij for neighboring 
atoms is negative. Such a magnetically ordered state is called antiferromagnetism.^EEHII!] 
In 1948, L. Neel introduced the notion of ferrimagnetisniH^HSS] ^ Q cl e g Cr ib e the properties of 
substances in which spontaneous magnetization appears below a certain critical temperature due 
to nonparallel alignment of the atomic magnetic moment liSSHUMSSl These substances differ from 
antiferromagnets where sublattice magnetizations and mg usually have identical absolute 
values, but opposite orientations. Therefore, the sublattice magnetizations compensate for each 
other and do not result in a macroscopically observable value for magnetization. In ferrimagnetics 
the magnetic atoms occupying the sites in sublattices A and B differ both in the type and in the 
number. Therefore, although the magnetizations in the sublattices A and B are antiparallel to 
each other, there exists a macroscopic overall spontaneous magnetization .HQSlI 1 1 1 II 1 12 ]H16JHI8J 
Later, substances possessing weak ferromagnetism were investigated .HOMES] n [ s interesting that 
originally Neel used the term parasitic ferromagnetisrrfr^ when referring to a small ferromagnetic 
moment, which was superimposed on a typical antiferromagnetic state of the a iron oxide Fe20$ 
(hematite) Later, this phenomenon was called canted antiferromagnetism, or weak ferromag- 
netismP^Hl^ 2 ^ The weak ferromagnetism appears due to antisymmetric interaction between the 
spins S\ and 52 and which is proportional to the vector product 5i x 52- This interaction is written 
in the following form 

H DM ~ DSi x 5 2 . (4.6) 

The interaction (|4,6p is called the Dzyaloshinsky-Moriya interaction [127, 128]. Hematite is one 
of the most well known minerals liMUlsglEfflMU which is still being intensively studiecP^ even at 
the present timeP^HESl 

Thus, there exist a large number of substances and materials that possess different types of 
magnetic behavior: diamagnetism, paramagnetism, ferromagnetism, antiferromagnetism, ferri- 
magnetism, and weak ferromagnetism. We would like to note that the variety of magnetism is 
not exhausted by the above types of magnetic behavior; the complete list of magnetism types is 
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substantially longerP^ As was already stressed, many aspects of this behavior can be reasonably 
well described in the framework of a very crude Heisenberg-Dirac-van Vleck model of localized 
spins. This model, however, admits various modifications (see, for instance, the boolf^). There- 
fore, various nontrivial generalizations of the localized spin models were studied. In particular, 
a modification of the Heisenberg model was investigated, where, in addition to the exchange in- 
teraction between different sites, an exchange interaction between the spins at the same site was 
considered?^ 

n = -vbh s l - \ E E J (*«; #o { xs t a sy p + sf a s^) - (4.7) 

<ia> i^j a/3 

\ E E J ^ #0 ( xs taSr P + sLsf?) ■ 

In the case when J(ia;if3) 3> J(ia;j{3), this model Hamiltonian in some sense imitates Hund's 
rules. Indeed, Hund's rules state that the triplet's spin state of two electrons occupying one and 
the same site is energetically more favorable than the singlet state. It is this feature that is taken 
into account by the model (|4,7p . A model of this type was used for description of composite 
ferrites, which contain different types of atoms with different spins (magnetic moments). In the 
limiting case J(ia;if3) ^ 0; J(i a; j(5 ) = the model (|4.7p can be considered as the simplest 
version of the Heisenberg modelfHSl i n this case, the two-spin system is interpreted^] a § the 
simplest one-dimensional periodic magnet with the period N = 2. Despite the apparent 'shortages' 
model (|4.7p has found numerous applications for description of real substances,^^ including the 
composite Cu(NOs)2 • 2.5i?20-type saltsp^ZEH Q f clusters ESSllHs] ag we rj ag f or improving mean- 
field approximation by using various cluster methodsP^Sl 

4.2 The Problem of Magnetism of Itinerant Electrons 

The Heisenberg model describing localized spins is mostly applicable to substances where the 
ground state's energy is separated from the energies of excited current-type states by a gap of 
a finite width. That is, the model is mostly applicable to semiconductors and dielectrics P-UJEZ] 
However, the main strongly magnetic substances, nickel, iron, and cobalt, are metals, belonging 
to the transition groupP^ The development of quantum statistical theory of transition metals 
and of their compounds followed a more difficult path than that of the theory of simple met- 
alsP^HHIH] The traditional physical picture of the metal state was based on the notion of Bloch 
electron waves D35HIS2 However, the role played by the inter-electron interaction remained un- 
clear within the conventional approach. On the other hand, the development of the band theory 
of magnetismj^^3Hl57| anc j investigations of the electronic phase's transitions in transition and 
rare-earth metal compounds gradually led to realization of the determining role of electron correla- 
tions PUSHES] Moreover, in many cases inter-electron interaction is very strong and the description 
in terms of the conventional band theory is no longer applicable. Special properties of transi- 
tion metals and of their alloys and compounds are largely determined by the dominant role of 
d-electrons. In contrast to simple metals, where one can apply the approximation of quasi-free 
electrons, the wave functions of (i-electrons are much more localized, and, as a rule, have to be 
described by the tight-binding approximationP^^^^ The main aim of the band theory of mag- 
netism and of related theories, describing phase ordering and phenomena of phase transition in 
complex compounds and oxides of transition and rare-earth metals, is to describe in the frame- 
work of a unified approach both the phenomena revealing the localized character of magnetically 
active electrons, and the phenomena where electrons behave as collectivized band entities.^ A 
resolution of this apparent contradiction requires a very deep understanding of the relationship 
between the localized and the band description of electron states in transition and rare-earth 
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metals, as well as in their alloys and compounds. The quantum statistical theory of systems with 
strong inter-electron correlations began to develop intensively when the main features of early 
semi-phenomenological theories were formulated in the language of simple model Hamiltonians. 
Both the Anderson model,' 162 " 163 ' which formalized the Friedel theory of impurity levels, and the 
Hubbard model JE2hl£2l which formalized and developed early theories by Stoner, Mott, and Slater, 
equally stress the role of inter-electron correlations. The Hubbard Hamiltonian and the Anderson 
Hamiltonian (which can be considered as the local version of the Hubbard Hamiltonian) play an 
important role in the electron's solid-state theory.^ Therefore, as was noticed by E. Liebj^the 
Hubbard model is 'definitely the first candidate' for constructing a 'more fundamental' quan- 
tum theory of magnetic phenomena than the 'theory based on the Ising model^™ (see also the 
paperimHIZS]). 

However, as it turned out, the study of Hamiltonians describing strongly-correlated systems is an 

exceptionally difficult many-particle problem, which requires applications of various mathematical 

methods EGllIZSES i n f ac t 5 with the exception of a few particular cases, even the ground state 

of the Hubbard model is still unknown. Calculation of the corresponding quasiparticle spectra in 

the case of strong inter-electron correlations also turned out to be quite a complicated problem. 

As was quite rightly pointed out by J. Kanamorip^ when one is dealing with "a metal state 

with the values of parameters close to the critical point, where the metal turns into a dielectric", 

then "the calculation of excited states in such crystals becomes very difficult (especially at low 

temperatures)". Therefore, in contrast to quantum many-body systems with weak interaction, 

the definition of such a notion as elementary excitations for strongly-interacting electrons with 

strong inter-electron correlations is quite a nontrivial problem requiring special detailed investiga- 
tionsPEJM] At the 

same time, one has to keep in mind, that the Anderson and Hubbard models 
were designed for applications to real systems, where both the case of strong and the case weak 
inter-electron correlations are realized. 

Often, a very important role is played by the electron interaction with the lattice vibrations, the 
phonons 11S2H1S2I Therefore, the number one necessity became the development of a systematic 
self-consistent theory of electron correlations applicable for a wide range of the parameter val- 
ues of the main model, and the development of the electron-phonon's interaction theory in the 
framework of a modified tight-binding approximation of strongly correlated electrons, as well as 
the examination of various limiting cases EMlM ^11 this activity allowed one to investigate the 
electric conductivity J^"'^ 6 -' and the superconductivity^^^ [ n transition metals, and in their 
disordered alloys. 

4.3 The Anderson and Hubbard Models 

The Hamiltonian of the single-impurity Anderson niodeP^lIIQSa] [ s written in the following form: 
ft = Yl e kc\ a c k(J + E 0a flf Qa + U/2 noano-a + ^ V k(4af0a + fl a c ka ), (4.8) 

ka a u kcr 

Here, c\ a and /q ct are the creation operators of conduction electrons and of localized impurity 
electrons, respectively, e k are the energies of conduction electrons, is the energy level of 
localized impurity electrons, and U is the intra-atomic Coulomb interaction of the impurity-site 
electrons; V k is the s — f hybridization. One can generalize the Hamiltonian of the single-impurity 
Anderson model to the periodic case: 

ft = E e kC ] ka C ka + E °fiJio + U / 2 E n - n *- + E V *3 (4a + (4-9) 

ka ia ia kja 
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The above Hamiltonian is called the periodic Anderson model. The Hamiltonian of the Hubbard 
modeP 31 is given by: 

H = J2 %4r°J> + U / 2 Yl n i<r n i-v ( 4 - 10 ) 
iju ia 

The above Hamiltonian includes the repulsion of the single-site intra-atomic Coulomb U, and tij, 
the one-electron hopping energy describing jumps from a J site to an i site. As a consequence of 
correlations electrons tend to "avoid one another". Their states are best modeled by atom-like 
Wannier wave functions [4>(f — Rj)]. The Hubbard model's Hamiltonian can be characterized by 
two main parameters: U, and the effective band width of tightly bound electrons 

A = (N~^\^\ 2 ) 1/2 - 

ij 

The band energy of Bloch electrons e(k) is given by 

e(k) = N' 1 ^2 Uj exp[-ik(Ri - Rj], 
k 

where N is the total number of lattice sites. Variations of the parameter 7 = A/U allow one to 
study two interesting limiting cases, the band regime (7 3> 1) and the atomic regime (7 — > 0). 
Note that the single-band Hubbard model (|4.10p is a particular case of a more general model, 
which takes into account the degeneracy of d-electrons. In this case the second quantization is 
constructed with the aid of the Wannier functions of the form [(j)\(r — Ri)], where A is the band 
index (A= 1,2,. . . 5). The corresponding Hamiltonian of the electron system is given by 

W = X) i ^ij a \ i ia a 3^ + 2 Y Yl (^JP\W\mj,nd)4 aa a t j p af a mitTl a nS(7 . (4.11) 

ijfj.ua ij,mn afj-yStTtj' 

It can be rewritten in the following form 

% = H l + H 2 + H z . (4.12) 
The first term here represents the kinetic energy of moving electrons 

ff i = EEC a U-. (4.13) 

ij fiver 

The second term H2 describes the single-center Coulomb interaction terms: 

ifia ifiv crcr' 



2 y / IfjLV n ifj,CF n ivcrO- — &fiv) + 2 ^ y ^ ^ a ificj a ifjL-cj a i^-cr a i^o{^ ~ ^ 

2 5-/ ^ V va ifia ai l J --v a iv—a a iva()- — <5^; 



Except for the integral of single-site repulsion f/«u, which is also present in the single-band Hubbard 
model, H2 also contains three additional contributions from the interorbital interaction. The last 
term H3 describes the direct inter-site's exchange interaction: 

ijfi a a' 
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It is usually assumed that 



U^ = U- V^ = V; V = J; J$* = Jy. (4.16) 

It is necessary to stress that the Hubbard model is most closely connected with the Pauli ex- 
clusion principle, which in this case can be written as nf a = rii a . Thus, the Anderson and the 
Hubbard models take into account both the collectivized (band) and the localized behavior of 
electrons. The problem of the relationship between the collectivized and the localized description 
of electrons in transition and rare-earth metals and in their compounds is closely connected with 
another fundamental problem. The case in point is the adequacy of the simple single-band Hub- 
bard model, which does not take into account the interaction responsible for Hund's rules and the 
orbital degeneracy for description of magnetic and some other properties of matter. Therefore, 
it is interesting to study various generalizations of the Anderson and the Hubbard models. In a 
series of paper^^^Sl we pointed out that the difference between these models is most clearly 
visible when we consider dynamic (as opposed to static) characteristics. Therefore, the response 
of the systems to the action of external fields and the spectra of excited quasiparticle states are 
of particular interest. Introduction of additional terms in the Anderson and the Hubbard model's 
Hamiltonians makes the quasiparticle spectrum much more complicated, leading to the appearance 
of new excitation branches, especially in the optical region ] 



4.4 The s-d Exchange Model and the Zener model 

A generalized spin-fermion model, which is also called the Zener model, or the s-d- (d-/)-model 
is of primary interest in the solid-state theory. The Hamiltonian of the s-d exchange modeP^ is 
given by: 

H = H S + H s _ d , (4.17) 
H s = ^e k c\ a c ka , (4.18) 

ka 

H s ^ d = JaiSi = -JN- 1 ' 2 (cUchS- + 4 a c fct S+ + (4, t c &t - c{, iCki )S z ) . (4.19) 

kk' 

Here, cj, and c kr7 are the second-quantized operators creating and annihilating conduction elec- 
trons. The Hamiltonian (14. 17ft describes the interaction of the localized spin of an impurity atom 
with a subsystem of the host-metal conduction's electrons. This model is used for description of 
the Kondo effect, which is related to the anomalous behavior of electric conductivity in metals 
containing a small amount of transition metal impurities ESlll2QHl92| 

It is rather interesting to consider a generalized spin-fermion model, which can be used for de- 
scription of a wider range of substances ESJUSHUSH] Hamiltonian of the generalized spin-fermion 
d-f modeP^l ig given by: 

H = H d + H d „ f , (4.20) 

H d = ^ ^ t ij a \a CL ^ + \ ^2 U n ™ n i-v ( 4 - 21 ) 
ij cr i(T 

The H d _f operator describes the interaction of a subsystem of strongly localized 4/(5/)-electrons 
with the spin density of collectivized d-electrons. 

H d-f = Yl JB ^ = ~ JN ~ 1/2 Yl a L a k+ q -<r + z a Sl q al a a k+qa }. (4.22) 

i kq a 
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The sign factors z a , introduced here for convenience, are given by 

= (+,-); -a = (t,|); S^ = \ S -?' ~ a = + ' (4.23) 

In the general case, the indirect exchange integral J depends significantly on the wave vector 
J(k;k + q) and attains the maximum value at the point k = q = 0. Note that the conduction 
electrons from the metal s-band are also taken into account by the model, and their role is 
the renormalization of the model parameters due to screening and other effects. Note that the 
Hamiltonian of the s-d model is a low-energy realization of the Anderson model. This can be 
demonstrated by applying the Schrieffer-Wolf canonical transformatioij ^ 192 !^^^ t G trie- latter 
model. 



4.5 Falicov-Kimball Model 

In 1969, Falicov and Kimball proposed a "simple" (in their opinion) model for description of 
the metal-insulator transition in rare-earth metal compounds. This model describes two subsys- 
tems, namely, the band and the localized electrons and their interaction with each other. The 
Hamiltonian of the Falicov-Kimball modeP^Sl [ s given by 

U = H + H inU (4.24) 

where 

^° = E E ^a+^tw + E E Eb l b ^'- ( 4 - 25 ) 

k vo i a 

Here, aJ, fe(J is the operator creating in the band v an electron in the state with the momentum k 
and the spin a, and b\ is the operator creating an electron (hole) with the spin a in the Wannier 
state at the lattice site R{. The energies e u {k) and E are positive and such that wm[E+e u (k)] > 0. 
It is assumed that due to screening effects only intra-atomic interactions play a significant role. 
Falicov and KimbalP^Sl took into account six different types of intra-atomic interactions, and 
described them by six different interaction integrals Gj. In a simplified mean-field approximation 
the Hamiltonian of the model (|4.24p was given by 

U = N[en a + En b - Gn a n b ], (4.26) 

where nj, = A r_1 Ylia ^LAo-- Then, one can calculate the free energy of the system, and to in- 
vestigate the transition of the first-order semiconductor-metal phase. The Falicov-Kimball model 
together with its various modifications and generalizations became very popular^2ZH2Q3j ^ n i nves ti_ 
gating various aspects of the theory of phase transitions, in particular, the metal-insulator transi- 
tion. It was also used in investigations of compounds with mixed valence, and as a crystallization 
model. Lately, the Falicov-Kimball model was used in investigations of electron ferroelectricity 
(EFE)PSl It also turned out that the behavior of a wide class of substances can be described 
in the framework of this model. This class includes, for instance, the compounds YblnCu^, 
EuNi2(Sii- x Ge x )2, NH2, Ta x N. Thus, the Falicov-Kimball model is a microscopic model of the 
metal-insulator phase's transition; it takes into account the dual band-atomic behavior of elec- 
trons. Despite the apparent simplicity, a systematic investigation of this model, as well as of the 
Hubbard model, is very difficult, and it is still intensively studied P^ZH2Q3] 
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4.6 The Adequacy of the Model Description 

As one can see, the Hamiltonians of s—d- and d—f- models especially, clearly demonstrate the man- 
ifestation of collectivized (band) and the localized behavior of electrons. The Anderson, Hubbard, 
Falicov-Kimball, and spin-fermion models are widely used for description of various properties 
of the transition and rare-earth metal compounds fHHSlEIlBSlEllMll i^ ll ^osrt inS] T n particular, they 
are applied for description of various phenomena in the chemical adsorption theorypSSI surface 
magnetism, in the theory of the quantum diffusion in solid He 3 , for description of vacancy motion 
in quantum crystals, and the properties of systems containing heavy fermions .IHlll221[l2Sll2iQii2iij 
The latter problem is especially interesting and it is still an unsolved problem of the physics 
of condensed matter. Therefore, development of a systematic theory of correlation effects, and 
description of the dynamics in the many-particle models (|4.8p - (|4.10p . (|4.17p and (|4.20p . were 
and still are very interesting problems. All these models are different description languages, 
different ways of describing similar many-particle systems. They all try to give an answer on the 
following questions: how the wave functions of, formerly, valence electrons change, and how large 
the effects of changes are; how strongly do they delocalize? Their applicability in concrete cases 
depends on the answers to those questions. On the whole, applications of the above mentioned 
models (and their combinations) allow one to describe a very wide range of phenomena and to 
obtain qualitative, and frequently quantitative, correct results. Sometimes (but not always) very 
difficult and labor-intensive computations of the electron band's structure add almost nothing 
essential to results obtained in the framework of the schematic and crude models described above. 
In investigations of concrete substances, transition and rare-earth metals and their compounds, 
actinides, uranium compounds, magnetic semiconductors, and perovskite-type manganites, most 
of the described above models (or their combinations) are used to a greater or lesser degree. This 
reflects the fact that the electron states, which are of interest to us, have a dual collectivized and 
localized character and can not be described in either an entirely collectivized or entirely localized 
form. As far back as 1960, C. HerringJ^EI in his paper on the d-electron states in transition metals, 
stressed the importance of a "cocktail" of different states. This is why efforts of many researchers 
are directed towards building synthetic models, which take into account the dual band-atomic 
nature of transition and rare-earth metals and their compounds. 

It was not by accident, that E. LietP^made the following statement: "Search for a model Hamil- 
tonian describing collectivized electrons, which, at the same time, is capable of describing correctly 
ferromagnetic properties, is one of the main current problems of statistical mechanics. Its impor- 
tance can be compared to such widely known recent achievements, as the proof of the existence 
of extensive free energy for macroscopically large systems" (see alsd^lEHll^lil) Solution of this 
problem is a part of the more general task of a unified quantum-statistical description of electrical, 
magnetic, and superconducting properties of transition and rare-earth metals, their alloys, and 
compounds. Indeed, the dual band-atomic character of d- and, to some extent, /-states manifests 
itself not only in various magnetic properties, but also in superconductivity, as well as in the 
electrical and thermal conduction processes. 

The Nobel Prize winner K.G. Wilson noticedPEl "There are a number of problems in science which 
have, as a common characteristic, a complex microscopic behavior that underlies macroscopic ef- 
fects" (see alscEEEHl). Eighty years since the formulation of the Heisenberg model (in 1928), 
we still do not have a complete and systematic theory, which would allow us to give an unam- 
biguous answer to the questionPISl " Why is iron magnetic 1 !" Although over the past decades the 
physics of magnetic phenomena became a very extensive domain of modern physics, and numer- 
ous complicated phenomena taking place in magnetically ordered substances found a satisfactory 
explanation, nevertheless recent investigations have shown that there are still many questions that 
remain without an answer. The model Hamiltonians described above were developed to provide 
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an understanding (although only a schematic one) of the main features of real-system behavior, 
which are of interest to us. It is also necessary to stress, that the two types of electronic states, the 
collectivized and the localized ones, do not contradict each other, but rather are complementary 
ways of quantum mechanical description of electron states in real transition and rare-earth metals 
and in their compounds. In some sense, all the Hamiltonians described above can be considered as 
a certain special extension of the Hubbard Hamiltonian that takes into account additional crystal 
subsystems and their mutual interaction. 

The variety of the available models reflects the diversity of magnetic, electrical, and superconduc- 
tivity properties of matter, which are of interest to us. We would like to stress that the creation 
of physical models is one of the essential features of modern theoretical physics According to 
Peierls, "various models serve absolutely different purposes and their nature changes accordingly 
. . . A common element of all these different types of models is the fact, that they help us to imag- 
ine more clearly the essence of physical phenomena via analysis of simplified situations, which are 
better suited for our intuition. These models serve as footsteps on the way to the rational expla- 
nation of real-world phenomena . . . We can take those models, turn them around, and most likely 
we would obtain a better idea on the form and structure of real objects, than directly from the 
objects themselves" The development of the physics of magnetic phenomena^^"^^2I] p roV es 
most convincingly the validity of Peierls' conclusion. 

5 Theory of Many-particle Systems with Interactions 

The research program, which later became known as the theory of many-particle systems with 
interaction, began to develop intensively at the end of 1950s - beginning of 60sPH Due to 
the efforts of numerous researchers: F. Bloch, H. Frohlich, J. Bardeen, N.N. Bogoliubov, H. 
Hugenholtz, L. Van Hove, D. Pines, K. Brueckner, R. Feynman, M. Gell-Mann, F. Dyson, R. 
Kubo, D. ter Haar, and many others, this theory achieved significant successes in solving many 
difficult problems of the physics of condensed matter J222H225| p^ ^ OQ ]p26] con t a ins an interesting 
details and the story about the development of some aspects of the theory of many-particle 
systems with interaction, and about its applications to solid-state physics. For a long time the 
perturbation theory (in its most diverse formulations) remained the main method for theoretical 
investigations of many-particle systems with interaction. In the framework of that theory, the 
complete Hamiltonian of a macroscopic system under investigation was represented as a sum of 
two parts, the Hamiltonian of a system of noninteracting particles and a weak perturbation: 

n = H + V. (5.1) 

In many practically important cases such approach was quite satisfactory and efficient. Theory of 
many particle systems found numerous applications to concrete problems, for instance, in solid- 
state physics, plasma, superfluid helium theory, to heavy nuclei, and many others. It is intensive 
development of the theory of many-particle systems that led to development of the microscopic 
superconductivity theory P23I22H1 Q u ite possibly, this was historically the first microscopic theory 
based on a sound mathematical foundationP22Hl32] development of the many-particle systems 
theory led to adaptation of many methods from quantum field theory to problems in statistical 
mechanics. Among the most important adaptations are the methods of Green's functions J222H22SI 
and the diagram technique. However, as the range of problems under investigation widened, 
the necessity to go beyond the framework of perturbation theory was felt more and more acutely. 
This became a pressing necessity with the beginning of theoretical investigations of transition and 
rare-earth metals and their compounds, metal-insulator transitions}^' and with the development 
of the quantum theory of magnetism. This necessity to go beyond the perturbation theory's 
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framework was felt by the founders of the Green's functions theory themselves. Back in 1951 J. 
Schwinger wroteP^S 

"... it is desirable to avoid founding the formal theory of the Green's functions on the 
restricted basis provided by the assumption of expandability in powers of coupling 
constants." 

Since the most important point of the theory of many-particle systems with interaction is an ade- 
quate and accurate treatment of the interaction, which can change (sometimes quite significantly) 
the character of the system behavior, in comparison to the case of noninteracting particles, the 
above remark by J. Schwinger seems to be quite farsighted. It is interesting to note, that, ap- 
parently, admitting the prominent role of J. Schwinger in development of the Green's functions 
method, N.N. Bogoliubov in his paper^^uses the term Green-Schwinger function (for an interest- 
ing analysis of the origin of the Green's functions method see the paper^Sl and also the boolP^). 
As far as the application of the Green's functions method to the problems of statistical physics is 
concerned, here, an essential progress was achieved after reformulation of the original method in 
the form of the two-time temperature Green's functions method. 



5.1 Two-time Temperature Green's Functions 

In statistical mechanics of quantum systems the advanced and retarded two-time temperature 
Green's functions (GF) were introduced by N.N. Bogoliubov and S.V. TyablikovP^S T n contrast 
to the causal GF, the above function can be analytically continued to the complex plane. Due to 
the convenient analytical property the two-time temperature GF is a very widespread method in 
statistical mechanics.^ 16 ^ " ' In order to find the retarded and advanced GF we have to use a 
hierarchy of coupled equations of motion together with the corresponding spectral representations. 
Let us consider a many-particle system with the Hamiltonian % = H — fiN; here p is the chemical 
potential and N is the operator of the total number of particles. If A(t) and B(t') are some 
operators relevant to the system under investigation, then their time evolution in the Heisenberg 
representation has the following form 



fmt\ ,,„ ( -/ 'Hi 



A(t) = exp I— j A(0) exp I— . (5.2) 

The corresponding two-time correlation function is defined as follows: 

(A(t)B(t')) = Tr(pA(t)B(t')), p = Z' 1 exp(-pH). 
This correlation function has the following property 

{A(t)B(t')) = (5.3) 
Z-'Tr (exp(- W >exp (f ) .4(0) exp (=^^) 5(0) exp (^)) 

= Z-'Tr (ex P (- W )exp ( m ~* ) ) -4(0) exp {^^) B(0)) 

= (A(t-t')B{0)) = (A{a)B(t' -t)). 

Usually it is more convenient to use the following compact notations (A(t)B) and {BA(t)), where 
t — t' is replaced by t. Since 

_ m + im = m+m (5 . 4) 
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these two correlation functions are related to each other. Indeed, we have 

(A(t)B) = (5.5) 
Z~ lr Tr (^exp(-/3^)exp (*™\ Aexp f 1 ^) exp(/3^) exp(-/3^)^ 

= Z~*» (ex P (- W)B exp ( M+M ) ^ eX p ( Z«g±M )) 

= (BA(t + ihf3)}. 

On can consider the correlation function (BA(t)) as the main one, because one can obtain the 
other function (A(t)B) by replacing the variable t — > t\ = t + ihf3 in (BA(t)). The spectral 
representation (Fourier transform over uj) of the function (BA(t)) is defined as follows: 

r+oo - 

(BA(t))= dujexp[--ut]J(B,A;u), (5.6) 

J-oo ri 

1 f +co i 
J(B,A;oj) = —J dtexp[-ut}(BA(t)). 

Equation (|5.6p is the spectral representation of the corresponding time correlation function. 
The quantities J(B,A;u) and J(A,B;u>) are the spectral densities (or the spectral intensities). 
It is convenient to assume that uj = huj c i as , where uj c ias is the classical angular frequency. The 
time correlation function can be written down in the following form 



(BA(t)) = Z- 1 ^{n\B\m){m\ exp[-Ht}Aexp[—Ht}\l)(l\ exp(-pH)\n) (5.7) 

nml 

= Z' 1 ^2(n\B\m) (m\A\n) exp(-/3e n ) exp ( ~j{^n - e m )t 



where 

U\n) = e n \n), exp[--%i]|rt) = exp(--e n t)\n) . 
Therefore, taking into account the identity 

1 f +co i 

-^^J dtexp[--(e n -e m -uj)t] = 5(e n -e m -uj), (5.8) 

we obtain 

J(B, A; uj) = Z- 1 ^2(n\B\m)(m\A\n) exp(-(3e n )5(e n - e m - uj). (5.9) 

nm 

Hence, the Fourier transform of the time correlation function is given by 

r+oo ; 

(A(t)B) = (AB(—t)) = / dujexp[-ujt]J(A,B;uj) (5.10) 

J-oo n 
r+oo i 

= / dujJ(A, B; — uj) exp[— — uit], 

J-oo ft 

where 

J (A, B; -uj) = Z~ l ^2(m\A\n)(n\B\m) exp(-/3e m )<5(e m - e n + uj) (5.11) 

nm 

= Z^ 1 ^2(n\B\m)(m\A\n) exp(—f3e n )5(e n — e m — uj) exp((3uj) 
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It is easy to check, that the following identity holds 

J(A, B; -co) = exp(0u)J(B, A; u). (5.12) 

For the spectral density of a higher order correlation function (B[A(t),7i\-) we obtain 

J(B,[A,H}-;u) =ojJ(B,A;u), (5.13) 
uJ(A,B;lo) = J(A, [H,B]_;oj) = J([A,H]-,B;w), 



Now we introduce the retarded, advanced, and causal GF: 

G r (A,B;t-t') = ((A(t),B(t'W = -ie(t-t')([A(t),B(t%),r, = ±, (5.14) 

G a (A, B; t-t') = ((A(t), B(t'))) a = i6{t' - t)([A(t), B(t%), r, = ±, (5.15) 

G C (A, B; t-t') = ((A(t), B{t'))) c = iT(A(t)B(t')} = (5.16) 
iO{t - t'){A(t)B(t')) + r/i6(t' - t)(B(t')A(t)), r) = ±. 

Here (. . .) is the average over the grand canonical ensemble, 9{t) is the Heaviside step function; 
the square brackets denote either commutator or anticommutator (r/ = ±): 

[A,B]-r, = AB-rjBA. (5.17) 

An important ingredient for GF application is their temporal evolution. In order to derive the 
corresponding evolution's equation, one has to differentiate GF over one of its arguments. Let us 
differentiate, for instance, over the first one, the time t. The differentiation yields the following 
equation of motion: 

id/dtG a {t, t') = 5(t - t')([A, B] v ) + (([A, H]-{t), B(t'))) a . (5.18) 

Here, the upper index a = r,a,c indicates the type of the GF: retarded, advanced, or causal, 
respectively. Because this differential equation contains the delta function in the inhomogeneous 
part, it is similar in its form and structure to the defining equation of Green's function from 
the differential equation theory^^ (about the George Green (1793-1841) creative activity see a 
detailed papeif^3). It is this similarity that allows one to use the term Green's function for the 
complicated object defined by Eqs. (|5.14p - (|5.17p . It is necessary to stress that the equations of 
motion for the three GF: retarded, advanced, and causal, have the same functional form. Only 
the temporal boundary conditions are different there. The characteristic feature of all equations 
of motion for GF is the presence of a higher order GF (relative to the original one) in the right 
hand side. In order to find the higher-order function, one has to write down the corresponding 
equation of motion for the GF (([A, H](t),B(t'))), which will contain a new GF of even higher 
order. Writing down consecutively the corresponding equations of motion we obtain the hierarchy 
of coupled equations of motion for GF. In principle, one can write down infinitely many of such 
equations of motion: 

(i) n d n /dt n G(t,t') = (5.19) 

n 

Y,d) n - k d n - k /dt n -H{t - f )<[[. • • [A,H]...H],B] V ) 
+(([[... [A,H]-...H}-(t),B(t'))). 
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The infinite hierarchy of coupled equations of motion for GF is an obvious consequence 
of interaction in many-particle systems. It reflects the fact that none of the particles (or, no 
group of interacting particles) can move independently of the remaining system. The next task 
is the solution of the differential equation of motion for GF. In order to do that one can use 
the temporal Fourier transform, as well as the corresponding boundary conditions, taking into 
account particular features of the problem under consideration. The spectral representation for 
GF, generalizing Eqs. (|5.6p - (|5.9p is given by 

/•oo ■ 

G r (A,B;t-t') = (27rh)- 1 dEG(A,B;E)exp[--E(t-t% (5.20) 

Joo 

POO ■ 

G(A,B;E) = ((A\B))e= dtG(A,B;t)exp(-Et). (5.21) 

Joo n 

On substitution of Eq. (|5.20|> in Eqs. f|5.18|) and (15.19R one obtains 

EG(A, B; E) = {[A,B] V ) + {{[A,H]-\B)) E ; (5.22) 

n 

E n G(A, B;E)=J2 ■ ■ [A W]. . .V], B] v ) (5.23) 

k=l 

+«[[... [A, n\- 

The above hierarchy of coupled equations of motion for GF (|5.23p is an extremely complicated 
and nontrivial object for investigations. Frequently it is convenient to rederive the same hierarchy 
of coupled equations of motion for GF starting from differentiation over the second time t' . The 
corresponding equations of motion analogous to Eqs. (|5.22p and (]5.23p are given by 

- EG(A, B; E) = -([A, B] v ) + ((A\ [B, (5.24) 

n 

(-l) n E n G(A, B-E) = - J2(-l) n ~ k E n ~ k ([A [...[B,n\... W]]„) 

k=1 k-l 

+((A\[...[B,H]-...H]-)) E . (5.25) 

n 

The main problem is how to find solutions of the hierarchy of coupled equations of motion for GF 
given by either Eq. (|5.23p or Eq. (|5.25|) ? In order to approach this difficult task one has to turn 
to the method of dispersion relations, which, as was shown in the papers by N.N. Bogoliubov and 
collaborators HI32I1243I j s q U ite an effective mathematical formalism. The method of retarded and 
advanced GF is closely connected with the dispersion relations technique,^ which allows one to 
write down the boundary conditions in the form of a spectral representation for GF. The spectral 
representations for correlation functions were used for the first time in the paper^El D y Callen and 
Welton (see alscP^) devoted to the fluctuation theory and the statistical mechanics of irreversible 
processes. GF are combinations of correlation functions 

f+oo ■ 

F A B(t - t') = (A(t)B(t')) = (A(t - t')B) = doj exp[-Lot}J(A, B; oo), (5.26) 

J — oo 



— oo 
oo 



r+oo ■ 

F B A(t' -t) = (B(t')A(t)) = (BA(t - t')) = du exp[--ut]J(B, A; u). (5.27) 

J — oo 

Therefore, the spectral representations for two-time temperature Green's functions can be written 
in the following form 



-oo 



£ — oj £ -co 
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where 

J'(B, A; u) = (exp(/3w) - rj)J(B, A; to) (5.29) 

and £ is the complex energy £ = Re£ + ilm£. 
Hence, 



\0J 



(5.30) 



/+oo / 
duo I J{B, A; u) exp(/3w) - tjJ(B, A; 
-oo \ 

= j + duo(j(B,A;—oo) —r)J(B,A;oj) \ = {AB — rjBA). 

Therefore, we obtain the following equation 

{{A\B)) £ = {AB - r/BA) + (([A,H]-\B)) e . (5.31) 

One should note that the two-time temperature Green's functions are not defined for t = t'; more- 
over, {{A(t)B(t'))) r = for t < t' ((A(t)B(t'))) a = t > t'. Using the following representations 
for the step- function 6(t) : 

9{t) = exp(-et)(e ->• 0, e > 0), t > 0; 6(t) = 0, t < 0. (5.32) 

we can rewrite the Fourier transform of the retarded (advanced) GF in the following form 

lim((A\B)} E±i£ = G r ^(A,B;E). (5.33) 

It is clear that the two functions, G r (A, B; E) and G a {A, B; E), are functions of a real variable 
E; they are defined as limiting values of the Green's function ((A\B))£ in the upper and lower 
half-plane, respectively. According to the Bogoliubov-Parasiuk theoremP^ISHSII] ^h e function 

= [ + °° ^AB.A^MM-v) {6U) 

J -oo £ — CO 

is an analytic function in the complex £-plane; this function coincides with G r (A, B; E) everywhere 
in the upper half-plane, and with G a (A, B; E) everywhere in the lower half-plane. It has singular- 
ities on the real axis; therefore, one has to make a cut along the real axis. Note that G r ^{A, B; t) 
is a generalized function in the Sobolev- Schwartz sense JHDI242H2MI f unc ti on G(A, B; £) is an 
analytic function in the complex plane with the cut along the real axis. It has two branches; one 
is defined in the upper half-plane, the other in the lower half-plane for complex values of £: 

mB))e J G '^ B ^ £ >°- (5.35) 
U ' " [G a (A,B;£), £<0. K ' 

The corresponding Fourier transform is given by 

/+oo • 
dEG r{a) (A,B;E)exp[-i:Et] = 
-oo « 

/+°° i r +OD T'(B A- 1, A 

„ ^-pi-i*' L ^tes- (5 - 36) 

Here, J' (B , A; oo) can be written down as follows (e —> 0) 

J'(B,A;oo) = ~ {{{A\B)) u+ie - {{A\B)) u . l£ ) . (5.37) 
2m 
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Therefore, the spectral representations for the retarded and the advanced GF are determined by 
the following relationships: 

G r (A,B;E) = ((A\B))l +l£ = (5.38) 

P % d \ ■ AB,A-u) 
J-oo E-uj + ie 

J-oo E — ui 

G a (A 1 B;E) = ((A\B))l_ i£ = (5.39) 

P % - ■ J\B,A;u) 
J-oo E-u-ie 

r+oo j'(b A;oj) 
= P dw v ' ' 1 +mJ'(B,A;E). 

J-oo E-U) 

In the derivation of the above equations we made use of the following relationshirj^^32H244| 

lim— !— -+P-^m8(x). (5.40) 

e^O X±IE X 

Here, P(l/x) indicates that one has to take the principal value when calculating integrals. As a 
result we obtain the following fundamental relationship for the spectral density 

_ lg(Ajj_GW) 
v 7 2vri exp(/?£) -77 v ; 

Thus, once we know the Green's function G r ^ (A, B; E) we can find J(B, A; E), and then calculate 
the corresponding correlation function. Using the relationship (I5.4ip . one can obtain the following 
dispersion relationships: 

^(W-^r^i^M. (5.42) 

71" J-OO E-Ul 

The most important practical consequence of the spectral representations for the retarded and 
advanced GF is the so called spectral theorem: 

(B(t')A(t)) = (5.43) 

1 f +OQ i 

/ dEexp[-E(t-t')][exp((3E)-7]}- 1 lmG A B(E + ie), 

t^J-oo n 

(A(t)B(t>)) = (5.44) 

1 f' + °° i 
/ dEexp(pE)exp[-E(t-t')][exp(pE) - ^ImG AB {E + ie). 

^ J-oo ft 

Equations (15.43H and (j5.44j) are of a fundamental importance for the entire method of two-time 
temperature GF. They allow one to establish a connection between statistical averages and the 
Fourier transforms of Green's functions, and are the basis for practical applications of the entire 
formalism for solutions of concrete problems .SDSimSEH 

5.2 The Method of Irreducible Green's Functions 

When working with infinite hierarchies of equations for GF the main problem is finding the meth- 
ods for their efficient decoupling, with the aim of obtaining a closed system of equations, which 
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determine the GF. A decoupling approximation must be chosen individually for every particular 
problem, taking into account its character. This "individual approach" is the source of critique 
for being too ad hoc, which sometimes appear in the papers using the causal GF and diagram 
technique. However, the ambiguities are also present in the diagram technique, when the choice 
of an appropriate approximation is made there. The decision, which diagrams one has to sum 
up, is obvious only for a narrow range of relatively simple problems. In the papers^^^ devoted 
to Bose-systems, and in the papers by the author of this review^2l'22H2l[26l[i93i254H256| ^gyoted to 
Fermi systems it was shown that for a wide range of problems in statistical mechanics and theory 
of condensed matter one can outline a fairly systematic recipe for constructing approximate so- 
lutions in the framework of irreducible Green's functions method. Within this approach one can 
look from a unified point of view at the main problems of fundamental characters arising in the 
method of two-time temperature GF. 

The method of irreducible Green's functions is a useful reformulation of the ordinary Bogoliubov- 
Tyablikov method of equations of motion. The constructive idea can be summarized as follows. 
During calculations of single-particle characteristics of the system (the spectrum of quasiparticle 
excitations, the density of states, and others) it is convenient to begin from writing down GF 
(15. 14ft as a formal solution of the Dyson equation. This will allow one to perform the neces- 
sary decoupling of many-particle correlation functions in the mass operator. This way one can 
to control the decoupling procedure conditionally, by analogy with the diagrammatic approach. 
The method of irreducible Green's functions is closely related to the Mori-Zwanzig's projection 
method J2SZH2S2I which essentially follows from Bogoliubov's idea about the reduced description of 
macroscopic systemsP^ In this approach the infinite hierarchy of coupled equations for correla- 
tion functions is reduced to a few relatively simple equations that effectively take into account the 
essential information on the system under consideration, which determine the special features of 
this concrete problem. 

It is necessary to stress that the structure of solutions obtained in the framework of irreducible 
GF method is very sensitive to the order of equations for GF^GH] j n which irreducible parts are 
separated. This in turn determines the character of the approximate solutions constructed on the 
basis of the exact representation. In order to clarify the above general description, let us consider 
the equations of motion ([5722]) for the retarded GF ([5714]) of the form {{A(t), A^(t'))) 

uG(u) = ([A,A%) + ({[A,B\-\A*)) U . (5.45) 

The irreducible (ir) GF is defined by 

W«L4,#]_|At)) = (([A,B\--zA\A% (5.46) 

The unknown constant z is found from the condition 

([^[A,H]_,A%) = 0. (5.47) 

In some sense the condition (j5.47j) corresponds to the orthogonality conditions within the Mori 
formalism J252H263] n [ s necessary to stress, that instead of finding the irreducible part of GF 
(( ir )(([A, H]_ | A* 1 ))) , one can absolutely equivalently consider the irreducible operators [A, £T]_) 
([A,H]^ ir \ Therefore, we will use both the notation ((A\B))) and (((A)^\B))), whichever 
is more convenient and compact. Equation (|5.47p implies 

{[A,A%) M„ 1 ' 

Here, Mq and Mi are the zero and first moments of the spectral density JS1|242JES1 Green's function 
is called irreducible (i.e. impossible to reduce to a desired, simpler, or smaller form or amount) 



24 



if it cannot be turned into a lower order GF via decoupling. The well-known objects in statistical 
physics are irreducible correlation functions (see, e.g. paperiP ^ 265 !). In the framework of the 
diagram techniqu e! 237 ^ the irreducible vertices are a set of graphs, which cannot be cut along 
a single line. The definition (]5.46p translates these notions to the language of retarded and 
advanced Green's functions. We attribute all the mean-field renormalizations that are separated 
by Eq. (15.46P to GF within a generalized mean field approximation 

G o = {[A,A\) (5 _ 49) 
[u-z) 

For calculating GF ([5^6]) . ^ (([A, H]_{t), A^ (t'))), we make use of differentiation over the second 
time t' . Analogously to Eq. (|5,46p we separate the irreducible part from the obtained equation 
and find 

G(u) = G°(oj) + GV)P((j)G°(u;). (5.50) 
Here, we introduced the scattering operator 

P = (M )- 1 (((([A^]-) (ir) |([^ t ,^]-) (ir) )))(Mo)- 1 . (5.51) 

In complete analogy with the diagram technique one can use the structure of Eq. (|5.50p to define 
the mass operator M: 

P = M + MG°P. (5.52) 

As a result we obtain the exact Dyson equation (we did not perform any decoupling yet) for 
two-time temperature GF: 

G = G° + G°MG. (5.53) 

According to Eq. (I5.52p . the mass operator M (also known as the self-energy operator) can be 
expressed in terms of the proper (called connected within the diagram technique) part of the 
many-particle irreducible GF. This operator describes inelastic scattering processes, which lead 
to damping and to additional renormalization of the frequency of self-consistent quasiparticle 
excitations. One has to note that there is quite a subtle distinction between the operators P 
and M. Both operators are solutions of two different integral equations given by Eqs. (I5.52p 
and (|5,53p . respectively. However, only the Dyson equation (|5.53p allows one to write down the 
following formal solution for the GF: 

G= [(G )- 1 -M]-\ (5.54) 

This fundamental relationship can be considered as an alternative form of the Dyson equation, 
and as the definition of the mass operator under the condition that the GF within the generalized 
mean- field approximation, G°, was appropriately defined using the equation 

G°G~ X + G°M = 1. (5.55) 

In contrast, the operator P does not satisfy Eq. (|5.55p . Instead we have 

(G )- 1 -G- 1 = PG°G-\ (5.56) 

Thus, it is the functional structure of Eq. (|5.54p that determines the essential differences 
between the operators P and M. To be absolutely precise, the definition (|5.52p has a symbolic 
character. It is assumed there that due to the similar structure of equations (j5.14j) - (I5.17P defining 
all three types of GF, one can use the causal GF at all stages of calculation, thus confirming the 
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sensibility of the definition (|5.52p . Therefore, one should rather use the phrase "an analogue of 
the Dyson equation". Below we will omit this stipulation, because it will not lead to misunder- 
standings. One has to stress that the above definition of irreducible parts of the GF (irreducible 
operators) is nothing but a general scheme. The specific way of introducing the irreducible parts 
of the GF depends on the concrete form of the operator A on the type of the Hamiltonian, and 
on the problem under investigation. 

Thus, we managed to reduce the derivation of the complete GF to calculation of the GF in the 
generalized mean-field approximation and with the generalized mass operator. The essential part 
of the above approach is that the approximate solutions are constructed not via decoupling of 
the equation-of-motion hierarchy, but via choosing the functional form of the mass operator in an 
appropriate self-consistent form. That is, by looking for approximations of the form M ~ 
Note that the exact functional structure of the one-particle GF (|5.54p is preserved in this approach, 
which is quite an essential advantage in comparison to the standard decoupling schemes. 

6 The Generalized Mean Fields 

Apparently, the mean field concept was originally formulated for many-particle systems (in an 
implicit form) in Van der Waals (1837-1923) Ph.D. thesis "On the Continuity of Gaseous and 
Liquid States". This classical paper was published in 1874 and became widely knownP*5l At 
first, Van der Waals expected that the volume correction to the equation of state would lead only 
to an obvious reduction of the available space for the molecular motion by an amount b equal 
to the overall volume of the molecules. However, the actual situation turned out to be much 
more complicated. It was necessary to take into account both corrections, the volume correction 
6, and the pressure correction a/V 2 , which led to the Van der Waals equation.^? Thus, Van 
der Waals realized that "the range of attractive forces contains many neighboring molecules". 
The development of this approach led to the insight, that one can try to describe the complex 
many-particle behavior of gases, liquids, and solids in terms of a single particle moving in an 
average (or effective) field created by all the other particles, considered as some homogeneous 
(or inhomogeneous) environment. That is, the many-particle behavior was reduced to effective 
(or renormalized) behavior of a single particle in a medium (or a field). Later, these ideas were 
extended to the physics of magnetic phenomena, where magnetic substances were considered 
as some kind of a peculiar liquid. That was the origin of the terminology magnetically soft 
and hard materials. Beginning from 1907 the Weiss molecular-field approximation^ became 
widespread in the theory of magnetic phenomena,^ and even at the present time it is still being 
used efficiently!^^ Nevertheless, back in 1965 it was noticed thatP^S 

"The Weiss molecular field theory plays an enigmatic role in the statistical mechanics 
of magnetism" . 

In order to explain the concept of the molecular field on the example of the Heisenberg ferromagnet 
one has to transform the original many-particle Hamiltonian (|4.4p into the following reduced one- 
particle Hamiltonian 

% = — 2/xo/^bS • h*" ^ . 
This transformation is achieved with the help of the identity 

S-S' = S-(S') + (S) ■ S' - (S) ■ (S') + C. 

Here, the constant C = (S — (S)) ■ (S' — (S')) describes spin correlations. The usual molecular- 
field approximation is equivalent to discarding the third term in the right hand side of the above 
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equation, and using the approximation C ~ (C) = (S ■ S') — (S) ■ (S 1 ). for the constant C. 



Let us consider this point in more detail. It is instructive to trace the evolution of the mean or 
concept of the molecular field for different systems. The list of some papers, which contributed to 
the development of the mean- field concept, is presented in Table 1. A brief look at that table allows 
one to notice a certain tendency. Earlier molecular-field concepts described the mean-field in terms 
of some functional of the average density of particles (n) (or, using the magnetic terminology, the 
average magnetization (M)), that is, as F[{n), (M)]. Using the modern language, one can say that 
the interaction between the atomic spins Si and their neighbors can be equivalently described by 
effective (or mean) field M m/ ). As a result one can write down 



Here, h ext is the external magnetic field, xo is the system's response function, and J(Rji) is the 
interaction between the spins. In other words, in the mean-field approximation a many-particle 
system is reduced to the situation, where the magnetic moment at any site aligns either parallel 
or anti-parallel to the overall magnetic field, which is the sum of the applied external field and 
the molecular field. 

Note that only the "averaged" interaction with i neighboring sites is taken into account, while the 
fluctuation effects are ignored. We see that the mean-field approximation provides only a rough 
description of the real situation and overestimates the interaction between particles. Attempts to 
improve the homogeneous mean-field approximation were undertaken along different directions 12ZQJ 
An extremely successful and quite nontrivial approach was developed by L. Neel J**^ who essentially 
formulated the concept of local mean fields (1932). Neel assumed that the sign of the mean-field 
could be both positive and negative. Moreover, he showed that below some critical temperature 
(the Neel temperature) the energetically most favorable arrangement of atomic magnetic moments 
is such, that there is an equal number of magnetic moments aligned against each other. This novel 
magnetic structure became known as the antiferromagnetism. 271 It was established that the 
antiferromagnetic interaction tends to align neighboring spins against each other. In the one- 
dimensional case this corresponds to an alternating structure, where an "up" spin is followed by a 
"down" spin, and vice versa. Later it was conjectured that the state made up from two inserted 
into each other sublattices is the ground state of the system (in the classical sense of this term). 
Moreover, the mean-field sign there alternates in the "chessboard" (staggered) order. 
The question of the true antiferromagnetic ground state is not completely clarified up to the 
present timeP^H22D This is related to the fact that, in contrast to ferromagnets, which have a 
unique ground state, antiferromagnets can have several different optimal states with the lowest 
energy. The Neel ground state is understood as a possible form of the system's wave function, 
describing the antiferromagnetic ordering of all spins.^^ Strictly speaking, the ground state is 
the thermodynamically equilibrium state of the system at zero temperature. Whether the Neel 
state is the ground state in this strict sense or not, is still unknown. It is clear though, that in the 
general case, the Neel state is not an eigenstate of the Heisenberg antiferromagnet's Hamiltonian. 
On the contrary, similar to any other possible quantum state, it is only some linear combination 
of the Hamiltonian eigenstates. Therefore, the main problem requiring a rigorous investigation is 
the question of Neel state's^ 7 - 7 -' stability. In some sense, only for infinitely large lattices, the Neel 
state becomes the eigenstate of the Hamiltonian and the ground state of the system. Nevertheless, 
the sublattice structure is observed in experiments on neutron scattering,^ and, despite certain 
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Table 1: The development of the mean-field concept 



Mean- field type 


Author 


Year 


A homogeneous molecular field in dense gases 


J.D. Van der Waals 


1873 


A homogenous quasi-magnetic mean-field 






in magnetics 


P.Weiss 


1907 


A mean- field in atoms: 






the Thomas-Fermi model 


L.H.Thomas, E.Fermi 


1926-28 


A homogeneous mean-field 






in many-electron atoms 


D.Hartree, V.A. Fock 


1928-32 


A molecular field in ferromagnets 


Ya. G. Dorfman, F.Bloch, 


1927 - 1930 


Inhomogeneous (local) mean-fields 






in antiferromagnets 


L.Neel 


1932 


A molecular field, taking into account 






the cavity reaction in polar substances 


L.Onsager 


1936 


The Stoner model of band magnetics 


E. Stoner 


1938 


Generalized mean-field approximation 






in many-particle systems 


T.Kinoshita, Y. Nambu 


1954 


The BCS-Bogoliubov mean-field 






in superconductors 


N.N. Bogoliubov 


1958 


The Tyablikov decoupling for ferromagnets 


S. V. Tyablikov 


1959 


The mean-field theory for the Anderson model 


P.W.Anderson 


1961 


The density functional theory for electron gas 


W.Kohn 


1964 


The Callen decoupling for ferromagnets 


H.B.Callen 


1963 


The alloy analogy (mean-field) 






for the Hubbard model 


J.Hubbard 


1964 


The generalized H-F approximation 






for the Heisenberg model 


Yu.A. Tserkovnikov, Yu.G. Rudoi 


1973-1975 


A generalized mean-field approximation 






for ferromagnets 


N.M. Plakida 


1973 


A generalized mean-field approximation 






for the Hubbard model 


A.L. Kuzemsky 


1973-2002 


A generalized mean-field approximation 






for antiferromagnets 


A.L. Kuzemsky, D. Marvakov 


1990 


A generalized random-phase approximation 






in the theory of ferromagnets 


A.Czachor, A.Holas 


1990 


A generalized mean-field approximation 






for band antiferromagnets 


A.L. Kuzemsky 


1999 


The Hartree-Fock-Bogoliubov mean-field 






in Fermi systems 


N.N. Bogoliubov, Jr. 


2000 
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worries,^ the actual existence of sublattices^^ is beyond doubt. 

Once Neel's investigations were published, the effective mean-field concept began to develop at a 
much faster pace. An important generalization and development of this concept was proposed in 
1936 by L. Onsager^^ in the context of the polar liquid theory. This approach is now called the 
Onsager reaction field approximation. It became widely known, in particular, in the physics of 
magnetic phenomenal 279 ^ 2 -^ In 1954, Kinoshita and NambrP^ developed a systematic method for 
description of many-particle systems in the framework of an approach which corresponds to the 
generalized mean-field concept. Later, various schemes of "effective mean-field theory taking 
into account correlations" were proposed (see the review^). One can show in the framework of 
the variation principle^^SlEHI] that various mean-field approximations can be described on the 
basis of the Bogoliubov inequality^ 

F = -/T 1 ln(Tre^) < 

-r'HTre-^ ) + >. (6.2) 

Here, F is the free energy of the system under consideration, whose calculation is extremely 
involved in the general case. The quantity H m f is some trial Hamiltonian describing the effective- 
field approximation. The inequality (|6.2p yields an upper bound for the free energy of a many- 
particle system. One should note that the BCS-Bogoliubov superconductivity theory222H230] - 1S 
formulated in terms of a trial (approximating) Hamiltonian, which is a quadratic form with respect 
to the second-quantized creation and annihilation operators, including the terms responsible for 
anomalous (or non-diagonal) averages. For the single-band Hubbard model the BCS-Bogoliubov 
functional of generalized mean fields can be written in the following forrrP^ 



r- I ( a i-a a i-°) -\<h««;-ni j - ., 

1 M-M ' ' ' 




The anomalous (or nondiagonal) mean values in this expression fix the vacuum state of the system 
exactly in the BCS-Bogoliubov form. A detailed analysis of Bogoliubov's approach to investiga- 
tions of (Hartree- Fock-Bogoliubov) mean-field type approximations for models with a four-fermion 
interaction is given in the papers E11286J 

There are many different approaches to construction of generalized mean-field approximations; 
however, all of them have a special-case character. The method of irreducible Green's functions 
allows one to tackle this problem in a more systematic fashion. In order to clarify this statement 
let us consider as an example two approaches for linearizing GF equations of motion. Namely, the 
Tyablikov approximation^ and the Callen approximation^^ for the isotropic Heisenberg model 
(|4.4p . We begin from the equations of motion ()5. 18[) for GF of the form {(S + \S~)): 

u>{{S?\Sr)) u = 2{S*)6ij g){(SfS x g - S+S*\Sr)) u . 

9 

Within the Tyablikov approximation the second order GF is written in terms of the first-order 
GF as followsP 

«#^|ST» * (sr)((s+\ST)). (6.4) 

It is well know, that the Tyablikov approximation (|6.4p corresponds to the random phase approx- 
imation for a gas of electrons. The spin-wave's excitation spectrum does not contain damping in 
this approximation: 

E(q) =J2 J ^~ 9)( S *) e*P[i(Ri ~ 4)3 = 2(^)(J - J q ). (6.5) 
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This is due to the fact that the Tyablikov approximation does not take into account the in- 
elastic quasiparticle's scattering processes. One should also mention that within the Tyablikov 
approximation the exact commutation relations [S^~,Sj]- = 2S z 5ij are replaced by approximate 
relationships of the form [S^,S~]- ~ 2(S z )5ij. Despite being simple, the Tyablikov approxima- 
tion is widely used in different problems even at the present timePSO 

Callen proposed a modified version of the Tyablikov approximation, which takes into account 
some correlation effects. The following linearization of equations-of- motion is used within the 
Callen approximationPS3 

((S* g S+\B)) (S z )((Sj\B)) - a(S;Sp((S+\B)). (6.6) 

Here, < a < 1. In order to better understand Callen's decoupling idea one has to take into 
account that the spin 1/2 operator S z can be represented in the form S z = S — SgS+ or S z = 
^(S^Sg — Sg S^). Therefore, we have 

1 — cv 1 -I- rv 

The operator S~Sg^ is the "deviation" of the quantity (S z ) from S. In the low-temperature 
domain that "deviation" is small and a ~ 1. Analogously, the operator ^(Sg^S~ — S~S£) is the 
"deviation" of the quantity {S z ) from 0. Therefore, when (S z ) approaches zero one can expect 
that a ~ 0. Thus, the Callen approximation has an interpolating character. Depending on the 
choice of the value for the parameter a, one can obtain both positive and negative corrections to 
the Tyablikov approximation, or even almost vanishing corrections. The particular case a = 
corresponds to the Tyablikov approximation. 

We would like to stress that the Callen approach is by no means rigorous. Moreover, it has 
serious drawbacks.^ However, one can consider this approximation as the first serious attempt 
to construct an approximating interpolation scheme in the framework of the GF's equations-of- 
motion method. In contrast to the Tyablikov approximation, the spectrum of spin-wave excitations 
within the Callen approximation is given by 

E(q) = 2(S*)((J - J q ) + ^J2i J ( k ) ~ J ( k ~ (6-7) 

k 

Here, N{E(k)) is the Bose's distribution function N(E(k)) = [exp(E(k)/3) - Equation (|57fjl 
clearly shows how the Callen approximation improves Tyablikov's approximation. From a general 
point of view, one has to find the form of the effective self-consistent generalized mean-field 
functional. That is, to find which averages determine that field 

F = {(S z ), (S x ), (S y ), (S + S~), (S Z S Z ), (s z s + s-), . . .}. 

Later many approximate schemes for decoupling the hierarchy of equations for GF were pro- 
posed,^ improving the Tyablikov and Callen decouplings. Various approaches generalizing the 
random phase's approximation in the ferromagnetism theory for wide ranges of temperature were 
considered in the paper^J by Czachor and Holas. 

6.1 Heisenberg Antiferromagnet and Anomalous Averages 

In order to illustrate the scheme of the irreducible GF method we are going to consider now the 
Heisenberg antiferromagnet. Note that a systematic microscopic theory of antiferromagnetism 
has not been built yet. In the framework of the model of localized spins the appearance of the 
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antiferromagnetic phase is usually associated with the first divergence of the generalized spin's 
susceptibility, if the exchange integral between the nearest neighbors is negative. The first di- 
vergence appears at Q = ir/a{a + b + c). Which means that when transiting from one atomic 
plane to another along the vector the phase of the magnetization vectors changes by ir. Generally 
speaking, in crystals with a complicated structure the exchange interaction may be different for 
different pairs of neighbors. In this case, we have a large variety of antiferromagnetic configura- 
tions. The simplest and the most frequently used model of localized spins of antiferromagnetic 
phenomena is the Heisenberg model of two-sublattice antiferromagnets. Let us consider now the 
calculation of the renormalized quasiparticle spectrum of magnetic excitations in the framework 
of the irreducible GF methodPSl The Hamiltonian of the system is given by 



% — — j)SiaSj a > — — SqaS-qa'- 

ij aa' q aa' 



(6.8) 



Here, Si a is the spin operator at the site i of the sublattice a, and J aa ' (i — j) is the exchange 
integral between the spins at the sites Ri a and Rj a >; the indexes a, a' assume two values (a) and 
(b). It is assumed that all the atoms in a sublattice a are identical and have the spin S a . It is 
convenient to rewrite the Hamiltonian (|6.8|) in the following form: 



H = - » EE ^ ( S qa S -qa' + S qa S -qa' ) > (6-9) 
q aa' 

where 

if = i/2(j« a ' + r' q a ). 

Let us again consider the equations of motion (|5.18j) for the Green's function of the form((5 + |S' _ )). 
In contrast to Heisenberg's ferromagnet model, for the two-sublattice antiferromagnet we have to 
use the matrix GF of the form 

GM - \((st b \sz ka )) ((st b \sz kb ))) • (6 - 10) 

Here, the GF on the main diagonal are the usual or normal GF, while the off-diagonal GF describe 
contributions from the so-called anomalous terms, analogous to the anomalous terms in the BCS- 
Bogoliubov superconductivity theory ()6.3|) . The anomalous (or off-diagonal) average values in this 
case select the vacuum state of the system precisely in the form of the two-sublattice Neel state. 
The Dyson equation (I5.53P is derived with the help of irreducible operators of the form 

(flg) <*•> = Sg - AfS+ a + A b k %S+, (6.11) 
(SfJ (-) = S z qa - N 1 ' 2 <S*> 6 q , , (6.12) 

where = (S^_ q a S qb — S qb Sf._ q a ). On performing standard transformations one can obtain the 
Dyson equation in the matrix form: 

G(k,u) = G (k,uj) + G (k,uj)M(k,uj)G(k,uj). (6.13) 

Here, Go(k,ui) is the GF within the generalized mean field approximation 

A _(G%*(k,u) Gg fe (fc,o;)\ _ 2<S z a > ((u-u aa ) oj ab \ 
Uo \G b a (k,co) G»(k,u,)J- det n V "a* (w-Wtt)J' 

where 

detCl = (U) - LO a a){u ~ ^bb) ~ ^aa^ab- 
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The poles of the GF (|6.14p determine the spectrum of magnetic excitations in the generalized 
mean- field approximation (the elastic scattering corrections): 



detQ = 0. 

As a result we obtain 



w±(fc) = ± y /(ul l (k)-u 2 ab (k)), (6.15) 

= [l - J^TjLl) E \/(l"7 fc 2 ), (6-16) 

where J 3 = 2/70, = 1/z exp(ifci?j) and z is the number of nearest neighbors. The first term 
in (|6.16p corresponds to the Tyablikov approximation. The second term describes the corrections 
of elastic scattering within the generalized mean-field approximation. Note that the quantity 
which determines these corrections, is given by 

. ab _ ns^ r \s^) + (sz ga s^ b ) 

A « - 2NW(S*) ' (6 ?) 

This expression contains anomalous averages (<S'Z 3a <S't ) , which characterize the Neel ground state. 

6.2 Many-particle Systems with Strong and Weak Electron Correlations 

The efficiency of the method of the irreducible Green's functions for description of normal and 
superconducting properties of systems with a strong interaction and complicated character of 
the electron spectrum was demonstrated in the papers l 2 ° ll 22 n 24 ll 254 l L e t, us consider the Hubbard 
model (14. lOH . The properties of this Hamiltonian are determined by the relationship between 
the two parameters: the effective band's width A and the electron's repulsion energy U. Drastic 
transformations of the metal-dielectric phase transition's type take place in the system as the ratio 
of these parameters changes. Note that, simultaneously, the character of the system description 
must change as well, that is, we always have to describe our system by the set of relevant variables. 
In the case of weak correlatiorP^J^HSlll^MI the corresponding set of relevant variables contains the 
ordinary second-quantized Fermi operators and a] a i(T) as well as the number of particles operator 
n ia = a \ a 0'i(7- These operators have the following properties: 

o t*(0) = tf(l). o .*(l)=«r(0) ) 

Here and vpW describe the vacuum and the single-particle states, respectively [159]. In order 
to find the low-lying excited quasiparticle states of the many-electron system with the Hamiltonian 
(|4.1U|) , one has to pass to the vector space of Bloch states 

= N~ 1/2 eyLp(-ikRi)a ia . 

i 

In this representation the Hamiltonian (|4.10p is given by 

(6.18) 

ka pqrs a 

Let us now consider the one-particle electron's GF of the form 

GW(i-0 = ((ak„,a{ a )) = -i9{t - t'){[a ka {t),a\ a (t')] + ) . (6.19) 
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The corresponding equation of motion (|5.18p for G ka (uj) is given by 

(u - e k )G ka {oj) = 1 + C/'/JV'5Z((afc+pff4+9-cr ?- <r l a L))w- ( 6 - 20 ) 



pq 

In line with Eq. ()5.46p we introduce the irreducible GF 

' ] {{ a k+pa( L \,J r q-a a q-(T\ a \a))^ = 

((ofc+paOp+g-^ a 9- CT l a L))^ ~ $p,o(n q -a)Gka- (6.21) 

The irreducible (ir) GF in Eq. (16,2ip is defined in such a way that it can not be transformed to a 
lower order GF by arbitrary pairings of second-quantized fermion operators. Next, according to 
Eqs. (prilD - (|534j) we find 

G ka (co) = GZ F (u) + (6.22) 
GZ F (u)U/N £ (*) ( {a k+pa al +q _ a a q ^ \ a\ a )) w . 

pq 

The following notation were introduced here 

G^ F (u) = (w - e(ka))- 1 ; e(ka) = e(k) + U/N^n^). (6.23) 

q 

Below, for simplicity we consider only paramagnetic solutions, where (n a ) = (n_ CT ). According to 
Eqs. (pT4~6|) - (f531j) we obtain 

G ka {u) = GH F (u) + G& F ( U )P^( U )G& F { U ). (6.24) 

The operator P is given by 

PkAuj) = ^ Y, D t ] (P> ^k, *,;«)= (6.25) 



^ { ( ir ) 1 1 t it t w(*»"A 

The proper part of the operator P is given by 

D ka(P^\r,s;u)) = L k "J(p,q\r,s;uj) 
+^ E LfJ{pMr'si^)Gti F {u)Dt\p\ q %s^). (6.26) 

r's'p'g' 

Here, L k %r J(p, q\r, s; oj) is the proper part of the GF D^(p,q\r,s;uj) Therefore, we obtain 

G ka = Gl'fiui) + Gff(u)M k(T (u)G K(T {u). (6.27) 

Equation (|6.27p is the desired Dyson equation for two-time temperature GF G ka (uj). It has the 
following formal solutions, cf. (|5.54p : 

G ka (u) = [u- e(ka) - M k(T {uS)]~ 1 . (6.28) 
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The mass operator M is given by 

M ktT (u) = = (6.29) 



U_ 

N 2 



pqrs 

2 / \ (P) 



■ W((a &+ p« T aJ +g _ CT a 9 _ CT |a;_ (T a r+s _ (T 4+ s <T)) (ir) J 



As was shown in the papers SOlEHMllSm an approximation to the mass operator M can be calculated 
as follows: 



P9 



X 



Cl> + Wl — UJ2 — W3 

n(w 2 )rt(o;3) + n(u)i)(l - n(w 2 ) - n(w 3 )^) J 5 p+9 _ (T (a;i)fi( A:+ p (T (a;2)5' (? - ( 7(w3). (6.30) 
Here, 

n(w) = [exp(/3w) + l] -1 ; # fcCT (u;) = lmG ka (u) + ie). 

Equations (|6.28p and (|6.30p are a self-consistent system of equations for calculating the one-particle 
GF Gfco-(w). As the first iteration one can substitute the expression 

9ka{u) ~ <$(w - e(fcer)). (6.31) 

in the right hand side of Eq. (|6.30p . The substitution yields 

i j / \ U 2 sr~^ np+g_ CT (l — rik+pa — riq- a ) + nk+ p(7 n q _ a 

M = ^ tr -+^+^)-^+^)-^) • ( } 

Equation (|6.32|) describes the renormalization of the electron spectrum due to the inelastic elec- 
tron's scattering processes. All elastic scattering corrections have already been taken into account 
by the electron energy's renormalization, see Eq. (|6.23p . Thus, the investigation of the Hubbard 
model in the weak coupling limit is relatively easy. 

The most challenging case is the solution of the Hubbard model when the electron correlations are 
strong, but are finite. In this limit it is convenient to consider the one-particle GF in the Wannier 
representations 

G ija (t - t') = ((a ta (t); a] a {t'))). (6.33) 

In the case of strong correlation, the algebra of relevant operators must be chosen according 
to specifi c fe atures of the problem under investigation. It is convenient to use the Hubbard 
operators?^ 

dicta — ^j_ CT fli(T, {pt — i), — f^iai ^ia — (1 ^i<r)j 

Yl n i° = X ; n i* n io = Sapn&i dicta = a ia . (6.34) 

a 

The new operators di au and dt„ have complicated commutation relations, namely 

[diaa,d j/3a ] + = 5 i j8 a i3nf_ u . 
The advantages of using these operators become clear when we consider their equations of motion: 

[dicta i H\— — Eadicta ~\~ ^ ^ tij(jl^_ (J Clj (7 + QMia-bij— CT ) , 

ij 

bija = {a\ a a ja - a^aia). (6.35) 
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According to Hubbardp^ the contributions to the above equation describe the "alloy analogy" 
corrections and the resonance broadening corrections. Using the Hubbard operators one can write 
down GF (|6.33p in the following form 

a/3 a/3 



The equation of motion for the auxiliary GF F 



is now given by 



F aP (u) 

11(7 \ / 



EF 8Jff (w) - 15, 



((di+a\4 +0 )) u ((d i+CT |4_ CT )) 6 



ij ) aft 



^tildni-o-atv + aa il7 b i i-. (r \d(p er )} ll 



Here, we used the following notation: 



E 



[u - E+) 






(oj - EJ 



n 







n 



(6.37) 



(6.38) 



(6.39) 



The determination of the irreducible parts of the GF is more involved: 



'((Zllldt-Hx))" {{Z 12 \d)_ a )) u 

S(Zn\d) +c )) u {{Z 22 \d)_ a )) u 



E 



A+a' 

a -a' 
■ A il 



\F a+ F a 

L ija iv 



Br 



i r lja r lja 



(6.40) 



In order to make the equations more compact we have introduced the following notation: 

Zu = Z12 = nf_ a ai a + aiabu^o-; Z 2 i = Z 22 = nl_ a ai a — ai a bii^ a . 

One has to stress that the definition (I6.40p plays the central role in this method. The coefficients 
A and B are found from the orthogonality condition (|5,47p 



{[(P { £3)aP,J jfh ,]+) 



0. 



(6.41) 



Next, the exact Dyson equation is derived according to Eqs. (|5.45p - (|5.54[) . Its mass operator is 
given by 



(6.42) 



I in 



The GF in the generalized mean-field's approximation has the following very complicated func- 
tional structurePl^HHEH] 



U) 



(n+E- + n_„E 4 



X(k) 



(6.43) 



(U3-E+- nI ff Ai(fc))(w -E-- nt a \ 2 {k)) - nZ a n± a \ 3 {k)\ A (k) 

Here, the quantities Aj(fc) are the components of the generalized mean field, which cannot be 
reduced to the functional of the mean particle's densities. The expression for GF (|6.43[) can be 
written down in the form of the following generalized two-pole solution 



n _ 



[l + cb- 1 ] 



a — db 1 c 



+ 



n 



(1 + da- 1 ) 



n 



u)-E. 



+ 



b — ca 1 d 
E + -n-_M-a 



(6.44) 



35 



where 



n± a nZ a W±_ a = N- 1 Uj exp[-ik(Ri - Rj)} x (6.45) 

ij 




Green's function (|6.44|) is the most general solution of the Hubbard model within the generalized 
mean field approximation. Equation (|6.45|) is nothing else but the explicit expression for the 
generalized mean field. As we see, this mean field is not a functional of the mean particle's 
densities. The solution (|6.44p is more general than the solution "Hubbard III"H2D anc i the two- 
pole solution from the papers^MSl D y R, th. It was shown in the papers^^HSlEHl by the author 
of this review, that the solution " Hubbard [ s a particular case of the solution (16. 44ft , which 
corresponds to the additional approximation 

n + anZaW ± {k) ra N'^tijeM-iKRi ~ Rj)]{nf^ntJ. (6.46) 

ij 

Assuming (nj~ a ni- a ) ~ n 2 _ a , we obtain the approximation "Hubbard I" .USD Thus, we have 
shown that in the cases of systems of strongly correlated particles with a complicated character of 
quasiparticle spectrums the generalized mean fields can have quite a nontrivial structure, which is 
difficult to establish by using any kind of independent considerations. The method of irreducible 
GF allows one to obtain this structure in the most general form. 



6.3 Superconductivity Equations 

The nontrivial structure of the generalized mean fields in many-particle systems is vividly revealed 
in the description of the superconductivity phenomenon. Let us now briefly consider this topic 
following the paperspSE2E3H223] "We describe our system by the following Hamiltonian: 

H = H e + H l + H e _i. (6.47) 

Here, the operator H e is the Hamiltonian of the crystal's electron subsystem, which we describe 
by the Hubbard Hamiltonian ()4.10p . The Hamiltonian of the ion subsystem and the operator 
describing the interaction of electrons with the lattice are given by 

^ = ^EH + ^E Ki«, (6-48) 

n mna/3 

= EE V%{Rl)a\ a a j(T u«, (6.49) 

where 

n ij 

Here, P n is the momentum operator, M is the ion mass, and u n is the ion displacement relative 
to its equilibrium position at the lattice site R n . Using more convenient notations one can write 
down the operator describing the interaction of electrons with the lattice as followd^^^^ 

H '-< = E E k + q)Q^a\ +qa a ka , (6.51) 

vo kq 



36 



where 

V u (k,k + q) = ^2K^a)e^(q)[sina a k- sina a (k- q)\. (6.52) 

Here, go is the Slater coefficient ESUEHIESH] describing the exponential decay of the d-electrons' wave 
function. The quantities e u (q) are the phonon- mode's polarization vectors. The Hamiltonian of 
the ion subsystem can be rewritten in the following form 

H * = \ E^^V + u 2 (qv)Q\ u Q qu ) (6.53) 

qv 

Here, P qv and Q qv are the normal coordinates, uj(qv) are the acoustic phonons' frequencies. 
Consider now the generalized one-electron GF of the following form: 

f r r \ I (( a ^l a jJ) ({aia\aj-a)) \ 

= (<£ I!) = m->i)) «aLi^» J • 

As was already discussed above, the off-diagonal entries of the above matrix select the vacuum 
state of the system in the BCS-Bogoliubov form, and they are responsible for the presence of 
anomalous averages. The corresponding equations of motion are given by 

y^juSjj - tij)((a ja \al a )) = 5u' + (6.54) 

3 

nj 

X>%+^)((a]->L>) = ( 6 - 55 ) 
3 

~U ({a\-<r n ia\4'*)) + v M{a)- a u n \a\, a )). 

nj 

Following the general scheme of the irreducible GF method, see Eqs. (|5.46p - (|5.54j) . we introduce 
the irreducible GF as follows 

{ {ir \(a ia a\_ a ai- a \a\ l(T ))u) = ((aurol_ ff a i _ .|4 (y )) u - (6.56) 
-(ni-a)Gxi + (a io .a i _ (r )((at_ ff |oT /o .)) w , 

( (<r) ((4^4->L>>«) = ({4* a iA-M>o))" - 

-{n i(T )G 2 l + (4a4-*)(( a ™\ a la))u- 

Therefore, instead of the algebra of the normal state's operator (a^, a\ a , rii a ) , for description of 
superconducting states, one has to use a more general algebra, which includes the operators (a,i a , 
olcr) n ia, a\<r a \-ai an d di-a<Ha)- The self-consistent system of superconductivity equations follows 
from the Dyson equation 

G iV {oj) = G!°,(a;) + J2 (^)Mjj' (w)Gj'i' (w) . (6.57) 

33' 

Green's function in the generalized mean- field's approximation, G°, and the mass operator Mj«i 
are defined as follows 

^2(u}T 5ij - tijT 3 - S^)G°j, = 6u>T , (6.58) 

3 
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M kk > = E(«(^ r 3^) M l(4 r 3P J 'fc') (jr) )))L p) , (6-59) 
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E 



M«/(w) = (6.60) 



^ V( (ir) ((«Wl^'n4't)) (ir) ) w ( (ir) ((4;^l^'^%-'^) (ir) ) (p V ' 

The mass operator ()6.60p describe the processes of inelastic electron scattering on lattice vibra- 
tions. The elastic processes are described by the quantity see Eq. (|6.3|) . An approximate 
expression for the mass operator (|6.60p follows from the following trial solution: 

{pj'i'a{t)a), a {t)a JCT p ija ) {vr) rj (pfi'cr{t)pi ja ) (4 ff (t)a j(T ). (6.61) 

This approximation corresponds to the standard approximation in the superconductivity theory, 
which in the diagram-technique language is known as neglecting vertex corrections, that is, ne- 
glecting electron correlations in the propagation of fluctuations of charge density. Taking into 
account this approximation, one can write down the mass operator (|6.60p in the following form 

M a iu) = Ml v (u>)+M&(u). (6.62) 

The first term^M 1 , has the form typical for an interacting electron-phonon system 

, ,i , . V^V^ T . rr 1 f + °° duj\duj2 ( 0U!l 0UJ2\ 

Mh h = e e 5 y ( cot ¥ + tan ¥ J x 

-~Im((it n |tv))w 2 J ^-^T 3 lm(('ipj\'i(j\ / )) ul T^j . (6.63) 
The second term M% has a more complicated structure 

Ml-, = % r d ^ ( cot ^ + tan , (6.64) 

2 y.^ w-wi-w2 v 2 2 / v m 2i ^ 22 y ' v ; 

where 



mn 

7T 



^Im((n 4 |n n )) W2 ^ f-^Im^a^lat^))^ J , (6.65) 
^Im((ra 4 |?v t )) W2 ^ ^-^Im^a^a^))^ , 
m 2i = ^Im((ni t |ni/|)) W2 ^ ^-^Im^a^a^))^ 
m 2 2 = y—lm^ni^rii^))^ f-^Im^a^lat^)), 



The definition (|6.56p and Eqs. (|6.57p - (|6.65p allowed us to perform a systematic derivation of 
superconductivity equations for transition met al ;p2l[23]lIMll293] anc j disordered binary alloys^^^El 
in the strong coupling approximation. Thus, it is the adequate description of the generalized 
mean-field in superconductors, taking into account anomalous mean values, which allowed us to 
construct compactly and self-consistently, the superconductivity equations in the strong coupling 
approximation. 
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6.4 Magnetic Polaron Theory 

To obtain a clear idea of the fundamental importance of the complex structure of mean fields 
let us investigate the problem of the magnetic polarorPSEHS - m ma g n etic semiconductors P33 
That is, in substances which have a subsystem of itinerant carriers and a subsystem of lo- 
cal magnetic moments EZll^MI|295] u sua lry the model of s — d exchange (|4.23p is used for de- 
scription of magnetic semiconductors. It is important to keep in mind that there are differ- 
ent spin and charge degrees of freedom in that model, which are described by the operators: 
a-ka, a\ a , n ka = a\ a a ka ; S£, SZ k = (S^; b ka = Y, q ( s ~q a q+k-u + z a S"L q a q+ka ); and 
a£ = Yl q a \^ a k+qi\ a k = Ylq a li a fc+<jt- The complete algebra of relevant operators is given 

Three additional GFs arise upon calculating the one-electron GF, because of the interaction be- 
tween the subsystems. In order to describe correctly the spin and charge degrees of freedom in 
magnetic semiconductors, as well as their interaction, the original GF must have the following 
matrix form: 



I ((a i<T \a^,)) «a»a|SJ» «<*to|S/'» «a j<7 |a^, S|» ((a i<T \a^_ a , sf » ^ 

((s?\J a ,)) «sfis|» «sf|sj'» ((sf\J gl sj)) «s;\<>l_ al sj')) 

((s-'w]^)) «s-"|s|» «s-"\sf)) ((s-"\a^,sp) «sr"l ]_„/ s i'» 

«Sfa i(T |a^,» «S?o i( ,|S|» «Sfa lo .|Sj'» «5f a ia \a* ol Sf)) «Sf *ia\*]_ a ,sf )) 

\((S-°a z _ a \a^,)) ((Si'ai-^Sf)) ((S-'a^^sf)) «Sr" a^^a^, Sf)) {{ sr^W]_ a ,sf))J 



(6.66) 



The functional structure of GF (|6.66|) shows that there are two regimes of quasiparticle dynamics: 
the scattering regime and the regime, where the electron-magnon's bound states (the magnetic 
polaron) are formed. To somewhat simplify our task we will use the following reduced algebra 
of relevant operators (a k(7 , a ka , b k(7 , b\). In this case, however, we will need a separate consistent 
consideration of the dynamic in the localized spin's subsystemPSlESSI For this purpose we use GF 

g + -(k;t-t') = ((S+(t),SZ k (t'))). (6.67) 

Now, the relevant matrix's GF for the problem of magnetic dynamics is given by 

mu) - {((4\sz k )) ((4\*-- k )))- (6 - 68) 

The Dyson equation for GF f)6.68j) 

g = So + g Mg (6.69) 

determines GF <3q in the generalized mean-field approximation, and the mass operator AfPSI For 
description of the charge-carriers subsystem we use the GF in the form 

9ka{t-t') = {{a ka (t),a{ a (t'))). (6.70) 

The Dyson equation for this GF is given by^S 

9ka(u) = g° ka (uj) + gl a {uj)M kcr (uj)g k(T (uj). (6.71) 

Equations (|6.69|) and (|6.7ip allow one to investigate self-consistently, the spin and the charge's 
quasiparticle dynamics in the system. In contrast to the scattering regime, for the one-electron 
GF (|6.70p in the bound state's formation regime we find the following expression for the GF in 
the generalized mean-field's approximation 

«a^|at»° = (detO)" 1 = (a, - e(ka) - ^Ar 1 *^))" 1 , (6.72) 
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where 

XM = (1 _ IAka{u)){uj + zMq) -e(k + q- a)) 

+ (l-/A feCT ( W ))(u;- e (fc + ga)) 1, 1 ' j 

N (oj + z a oj(q) - e(k + q- a)) 

The quantity plays the role of the generalized susceptibility for spin-electron bound states. 

It is this property that distinguishes the bound-state regime from the scattering regime, where 
instead of the electron-spin susceptibility x\ a ( u ) appears Xo(k,oj) 

X s (k,u) = iV-^ fe^zM , (6.75) 

p U P,k 

We use the following notation 

^p,k = (uj + e p - e p+k - A/); A/ = 2IS Z , 

i q <? 

e(qa) =e q - z a IS z , 



n = ^(n t + n ; ); < n < 2; 5, = AT 1 / 2 ^*). 

. The magnetic polaron's spectrum is given by 

E fc(y = e(fcff) + I 2 N- l X b ka (E ka ). (6.76) 

One can show that for any value of the electron's spin projection the polaron spectrum of the 
bound electron-magnon's state contains two branches. In the so-called atomic limit (e k = ), 
when k — > 0, oj — > 0, we obtain 

«a*>L»° = ^sTT^ + is)- 1 + ^sTT^ ~ 1{s + (6 ' 77) 

Here, S and S z = (<Sq) /y/~N denote the spin magnitude and the magnetization, respectively. The 
obtained result, Eq. (|6.76p . is in perfect agreement with the result of Mattis and ShastrypSS who 
investigated the magnetic polaron's problem for T = 

((a ka \al)) \ T =o = W- e{ka) - 5 ai 2I 2 S A *; M j' 1 (6.78) 

(1 - lA kcr (oj)) 

Thus, the magnetic polaron is formed in the case of antiferromagnetic s-d interaction (7 < 0). In 
order to get a clear idea of the spectrum character let us now consider two limiting cases: 

(i) a wide-band semiconductor (\I\S <C W) 

T S(S + S Z + 1) + S Z (S-S Z + 1) , . 

E K ^e k + 7^ - '-^^ - '- + (6.79) 



7) ^ (e k . q -e k + 2I(S - S z )) (S+S_ q ) 



N ^ (e k - q -e k + 2IS z ) 2S 
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(ii) a narrow-band semiconductor (\I\S 3> W) 

E H ~ 7(5 + 1) + (25 + 1)(5 + 5a + 1) e* + 

1 ( £ fc-9 ~ £ fc) ( S q S -q) ffi8m 

TV ^ (25 + 1) (5 + 5* + 1) ' { ' } 

Here, W is the band width for 1 = 0. Note that in order to make expressions more compact we 
omitted the correlation function K* z in the above formulae. 

Consider now the low-temperature spin-wave regime, where one can assume that S z ~ 5. In this 
case we have 

(S+SZ q )^2S(l+N(u(q))). 

One can show that for 

(i) a wide-band semiconductor (\I\S <C W) 

Ta 27 2 5^ 1 

N ^ (e k - e k _ q + 275) 

(ii) a narrow-band semiconductor (|7|5 S> W) 

Eki a I(S + 1) + ^ )tk + 1 2 ^j^Zgl ^W). ( 6 . 82 ) 

Let us now estimate the energy of the bound state of the magnetic polaron 

e B = e H - E H . (6.83) 

Taking into account that 

£kl = £fe + 75. 

we obtain the following expressions for the binding energy eb- 

(i) a wide-band semiconductor (|7|5 <C W) 

'' = ^-^E ( JX-*V ( " W) - (6 - 84) 

(ii) a narrow-band semiconductor (|7|5 S> W) 
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_ _ (2J 2 S) V 1 |I|5 

' B1 2V Y (e fc _ 9 - e fc - 2/5) W 1 

(25 + 1) 



B2 = -^+7^r L ^-m- (6-86) 



The outlined theory gives a complete description of the magnetic polaron for finite tempera- 
£ ureg |294][295] revea li n g the fundamental importance of the complicated structure of generalized 
mean-fields, which cannot be reduced to simple functionals of mean spin and particle densities. 



7 Broken Symmetry, Quasiaverages, and Physics of Magnetic 
Materials 

It is well known that the concept of spontaneously broken symmetry^ZHSS] j s one Q f the most 
important notions in the quantum field theory and elementary particle physics. This is especially 
so as far as creating a unified field theory, uniting all the different forces of naturep^ is concerned. 
One should stress that the notion of spontaneously broken symmetry came to the quantum field 
theory from solid-state physics. It was originated in quantum theory of magnetism, and later was 
substantially developed and found wide applications in the gauge theory of elementary particle 
physicsPMII] It was in the quantum field theory where the ideas related to that concept were 
quite substantially developed and generalized. The analogy between the Higgs mechanism giving 
mass to elementary particles and the Meissner effect in the Ginzburg-Landau superconductivity 
theory is well known p92ll298j[30i}t303Ji306Jl3ioj B th effects are consequences of spontaneously broken 
symmetry in a system containing two interacting subsystems. 

A similar situation is encountered in the quantum solid-state theory Analogies between the 
elementary particle and the solid-state theories have both cognitive and practical importance for 
their development i 3 -^ We have already mentioned the analogies with the Higgs effect playing an 
important role in these theories J 3 -^ However, we have every reason to also consider analogies with 
the Meissner effect in the Ginzburg-Landau superconductivity model, because the Higgs model 
is, in fact, only a relativistic analogue of that model 30 1 tt 3Q3J1306j Q n ^ e same ground one can 
consider the existence of magnons in spin systems at low temperatures}^] acoustic and optical 
vibration modes in regular lattices or in multi-sublattice magnets, as well as the vibration spectra 
of interacting electron and nuclear spins in magnetically-ordered crystals! 3 -^ 
The isotropic Heisenberg ferromagnet (|4.4|) is often used as an example of a system with spon- 
taneously broken symmetry! 3 -^ This means that the Hamiltonian symmetry, the invariance with 
respect to rotations, is no longer the symmetry of the equilibrium-state. Indeed the ferromag- 
netic states of the model are characterized by an axis of the preferred spin alignment, and, hence, 
they have a lower symmetry than the Hamiltonian itself. However, as was stressed by Ander- 
gon |3iiJ[3i6l|3i7| the ground state of the Heisenberg ferromagnet is an eigenstate of the relevant 
transformation of continuous symmetry (spin rotation). Therefore, the symmetry is not bro- 
ken and the low-energy excitations do not have novel properties. The symmetry breaking takes 
place when the ground state is no longer an eigenstate of a particular symmetry group, as in anti- 
ferromagnets or in superconductors. Only in this case the concepts of quasi-degeneracy, Goldstone 
bosons, and Higgs phenomenon can be applied EDlMMSI] 

The essential role of the physics of magnetism in the development of symmetry ideas was noted in 
the paper^^by the 2008 Nobel Prize Winner Y. Nambu, devoted to the development of the elemen- 
tary particle physics and the origin of the concept of spontaneous symmetry breakdown. Nambu 
points out that back at the end of the 19th century P. Curi c 319 ^ 32 ™ used symmetry principles in 
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the physics of condensed matter. P. Curie^^ used symmetry ideas in order to obtain analogues of 
selection rules for various physical effects, for instance, for the Wiedemann effedP^22Q] ( see the 
books^SHES^ Nambu also notes: 

"More relevant examples for us, however, came after Curie. The ferromagnetism is the prototype 
of today's spontaneous symmetry breaking, as was explained by the works of Weiss Heisen- 
bergGSD anc j others. Ferromagnetism has since served us as a standard mathematical model of 
spontaneous symmetry breaking" . 

This statement by Nambu should be understood in light of the clarification made by Ander- 
SOI] |3iij|3i6]|3i7| ( see a i so the papeJ^J). P. Curie was indeed a forerunner of the modern concepts of 
the quantum theory of magnetism. He formulated the Curie principle: "Dissymmetry creates 
the phenomenon" . According to this principle j2!9Jl320j 

"A phenomenon can exist in a medium possessing a characteristic symmetry (Gi) or 
the symmetry of one of that characteristic symmetry subgroups (G C Gi)" . 

In other words, some symmetry elements may coexist with some phenomena, but this is not 
necessarily the case. What is required is that some symmetry elements are absent. This is 
that dissymmetry, which creates the phenomenon. One of the formulations of the dissymmetry 
principle has the following forrrl 3 - 2 ^ 

^phenomena -j q ^phenomena 

or, alternatively 

^properties s~t f \ ^properties (r* n \ 

Lrj =i ^object — \ ■ \' 

Note that the concepts of symmetry, dissymmetry, and broken symmetry became very widespread 
in various branches of science and artJ^JEUD 

Essential progress in the understanding of the spontaneously broken symmetry concept is con- 
nected with Bogoliubov's ideas about quasiaveragesJ^^EMl 

Indeed, as was noticed in the book! 3 -^ _ .the canonical ensemble p ~ exp(—/3H) is no longer a 
good ensemble for the spontaneously ordered systems. Averaging over this ensemble would be to 
average, among other properties, over all directions of the total spin. That is fine in a paramagnet, 
and passes for a number of purposes in the ferromagnetic regime as well, but for other purposes, 
such as the calculation of {S tot }, it would be a foolish thing to do. One could use exp(—(3H) to 
weight states of different energy, but in addition one should specify that the trace is to be taken 
only over those states for which S tot points in the z-direction. Formally, one would then have 
something like 

p = const VS tot exp(-PH), 

where the projection operator VS tot eliminates all but those states for which S tot points along z." 
As we see, this statement written in 1975 contains in a concise form an argumentation in favor 
of using the ideas of quasiaver ages J323l328| ^ does not mention them explicitly. However, the 
notion of quasiaverage d 32 ' H 328 l was formulated by N.N. Bogoliubov back in 1960-1961 (see also the 
papeP31). 

It is necessary to stress, that the starting point for Bogoliubov's paper^^ 3 -^ was an investigation 
of additive conservation laws and selection rules, continuing and developing the already 
mentioned above approach by P. Curie for derivation of selection rules for physical effects. Bo- 
goliubov demonstrated that in the cases when the state of statistical equilibrium is degenerate, as 
in the case of a ferromagnet, one can remove the degeneracy of equilibrium states with respect to 
the group of spin rotations by including in the Hamiltonian H an additional noninvariant term 
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vM z V with an infinitely small v. This replaces the ordinary averages by quasiaverageJ^3G2S of 
the form 

lA^=1im(A) u g, (7.3) 

v — ^0 

where {A) u g is the ordinary average of the quantity A with respect to the Hamiltonian H u g = 
H + v{e- M)V. Thus, the presence of degeneracy is directly reflected on quasiaverages via their de- 
pendence on the arbitrary vector e. The ordinary averages can be obtained from the quasiaverages 
by integrating over all possible directions of e 

{A) = J <A> de. (7.4) 

The question of symmetry breaking within the localized and band models of antiferromagnets 
was studied by the author of this review in the paperspQEISlGSQ] j^. been found there that the 
concept of symmetry breaking in the band model of magnetisirPSS is much more complicated than 
in the localized model. In the framework of the band model of magnetism one has to additionally 
consider the so-called anomalous propagators of the form 

FM:G /m ~ ((a fcCT ;4_ CT )), 
AFM : G afm ~ ((a k+Qa ;al +Q/a ,}}. 

In the case of the band antiferromagnet the ground state of the system corresponds to a spin- 
density wave (SDW), where a particle scattered on the internal inhomogeneous periodic field gains 
the momentum Q — Q' and changes its spin: a — > a'. The long-range order parameters are defined 
as follows 



FM: m = l/Nj2(4a a k- C r}, (7.5) 
AFM : M Q = J2(4a a k+Q-a). (7.6) 

It is important to stress, that the long-range order parameters here are functionals of the internal 
field, which in turn is a function of the order parameter. Thus, in the cases of rotation and trans- 
lation invariant Hamiltonians of band ferro- and antiferromagnetics one has to add the following 
infinitesimal sources removing the degeneracy: 

FM : vn B H x a\a a k-v, (7-7) 

ka 

AFM : ufi B H ^ a\ a a k+Q ^. (7.8) 

kQ 

Here, v — > after the usual in statistical mechanics infinite-volume limit V —> oo. The ground state 
in the form of a spin-density wave was obtained for the first time by Overhauser in investigations 
of nuclear matterF^S] There, the vector Q is a measure of inhomogeneity or translation symmetry 
breaking in the system. It was written in the papei®2l ( see a i sc pSIH333^ 

"... in antiferromagnets a staggered magnetic field plays the role of a symmetry- 
breaking field. No mechanism can generate a real staggered magnetic field in antifer- 
romagnetic materials" . 

The analysis performed in the papers by p enr ]224||335j s h wed (see also^EI) that the antiferromag- 
netic and more complicated states (for instance, ferrimagnetic) can be described in the framework 



44 



of a generalized mean-field approximation. In doing that we have to take into account both the 
normal averages {a\ a a,ia)i and the anomalous averages (a/ ia ai- a ). 

It is clear that the anomalous terms (|7.7p and f|T.8j) break the original rotational symmetry of 
the Hubbard Hamiltonian. Thus, the generalized mean-field's approximation has the following 
form rij-o-ajo- ~ (nj_ -)aj (7 — (a\_ a ai a )ai- a . A self-consistent theory of band antiferromagnetism 
was developed by the author of this review in the paperP^^ using the method of the irreducible 
GF. The following definition was used: 

{{ a k+pa a p+q-a a q-o-\ a \ :cr ))uJ = {{ a k+pcr a \+q- (J a q-a\ a \ a )) u ~ 

$ P ,o{n q - a )G ka - (a fc+p(7 aJ +g _ (7 )((ag_ (T |4 (T )>a;- (7.9) 

The algebra of relevant operators must be chosen as follows (ai a , aj CT , rii a , a\ a ai- a ). The corre- 
sponding initial GF will have the following matrix structure 

qAFM ^ ( (( a iv\ a ]-a)) \ 

\((ai-a\a] a }} {{a^ a \a)_ a )) ) 

The off-diagonal terms select the vacuum state of the band's antiferromagnet in the form of a spin- 
density wave. It is necessary to stress that the problem of the band's antiferromagnetisuP^^ 
is quite involved, and the construction of a consistent microscopic theory of this phenomenon 
remains a topical problem. 



7.1 Quantum Protectorate and Microscopic Models of Magnetism 

The "quantum protectorate" concept was formulated in the paper P2D j^g authors, R. Laughlin 
and D. Pines, discussed the most fundamental principles of matter description in the widest sense 
of this word: 

"It is possible to perform approximate calculations for larger systems, and it is through such 
calculations that we have learned why atoms have the size they do, why chemical bonds have the 
length and strength they do, why solid matter has the elastic properties it does, why some things 
are transparent while others reflect or absorb light. With a little more experimental input for 
guidance it is even possible to predict atomic conformations of small molecules, simple chemical 
reaction rates, structural phase transitions, ferromagnetism, and sometimes even superconducting 
transition temperatures. But the schemes for approximating are not first-principles deductions 
but are rather art keyed to experiment, and thus tend to be the least reliable precisely when reli- 
ability is most needed, i.e., when experimental information is scarce, the physical behavior has no 
precedent, and the key questions have not yet been identified. . . .We have succeeded in reducing 
all of ordinary physical behavior to a simple, correct Theory of Everything only to discover that 
it has revealed exactly nothing about many things of great importance.' 

R. Laughlin and D. Pines show that there are facts that are clearly true, (for instance, the value 
e2/hc) yet they cannot be deduced by direct calculation from the Theory of Everything, for exact 
results cannot be predicted by approximate calculations. Thus, the existence of these effects is 
profoundly important, for it shows us that for at least some fundamental things in nature the The- 
ory of Everything is irrelevant. Next, the authors formulate their main thesis: emergent physical 
phenomena, which are regulated by higher physical principles, have a certain property, typical for 
these phenomena only. This property is their insensitivity to microscopic description. Thus, here, 
in essence, a most broad question is posed: 
"what is knowable in the deepest sense of the term?" 
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For instance, the low-energy excitation spectrum of ordinary crystal dielectrics contains a 
transversal and longitudinal sound wave and nothing else, irrespective of microscopic details (see 
alscP^). Therefore, in the opinion of R. Laughlin and D. Pines, there is no need "to prove" the 
existence of sound in solid bodies; this is a consequence of the existence of elastic modules in the 
long-wave scale, which in turn follows from the spontaneous breaking of translation and rotation 
symmetries, typical for the crystal state. This implies the converse statement: very little one 
can learn about the atomic structure of the solid bodies of crystal by investigating their acoustic 
properties. Therefore, the authors summarize, the crystal state is the simplest known example of 
the quantum protectorate, a stable state of matter with low-energy properties deter- 
mined by higher physical principles and by nothing else. 

The existence of two scales, the low-energy and high-energy scales, relevant to the description of 
magnetic phenomena was stressed by the author of this review in the papers^^'^ 8 - 9 -' devoted to 
comparative analysis of models of localized and band models of quantum theory of magnetism. It 
was shown there, that the low-energy spectrum of magnetic excitations in the magnetically-ordered 
solid bodies corresponds to a hydrodynamic pole (k,uj — > 0) in the generalized spin susceptibility, 
which is present in the Heisenberg, Hubbard, and the combined s — d model (see Fig. 1). 
In the Stoner band model the hydrodynamic pole is absent, there are no spin waves there. At the 
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Figure 1: Schematic diagrams of excitation spectra in four microscopic models of theory of mag- 
netism. Upper left: the Heisenberg model; upper right: the Hubbard model; lower left: the Zener 
model; lower right: the multiband Hubbard model. 

same time, the Stoner single-particle's excitations are absent in the Heisenberg model's spectrum. 
The Hubbard mode l 18 ^ 19 ^ 189 ! with narrow energy bands contains both types of excitations: the 
collective spin waves (the low energy spectrum) and Stoner single-particle's excitations (the high- 
energy spectrum). This is a big advantage and flexibility of the Hubbard model in comparison 
to the Heisenberg model. The latter, nevertheless, is a very good approximation to the realistic 
behavior in the domain where the hydrodynamic description is applicable, that is, for long wave- 
lengths and low energies. The quantum protectorate conc ept was applied to the quantum theory 
of magnetism by the author of this review in the paperf^SI where a criterion of models of the 
quantum theory of magnetism applicability to description of concrete substances was formulated. 
The criterion is based on the analysis of the model's low energy and high-energy spectra. 
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8 The Lawrence-Doniach Model 



The Ginzburg-Landau model 3 -^" 3 - 3 -^ is a special form of the mean-field theory. This model operates 
with a pseudo-wave function \l/(r), which plays the role of a parameter of complex order, while 
the square of this function modulus ^(r)! 2 should describe the local density of superconducting 
electrons. It is well known, that the Ginzburg-Landau theory is applicable if the temperature of 
the system is sufficiently close to its critical value T c , and if the spatial variations of the functions 
\P and of the vector potential A are not too large. The main assumption of the Ginzburg-Landau 
approach is the possibility to expand the free-energy density / in a series under the condition, 
that the values of are small, and its spatial variations are sufficiently slow. Then, we have 



/ = /, 



1 



2m* 



ihV H I * 



h 2 



The Ginzburg-Landau equations follow from an applications of the variational method to the 
proposed expansion of the free energy density in powers of |^| 2 and |V^| 2 , which leads to a pair 
of coupled differential equations for \P(r) and the vector potential A. 

The Lawrence-Doniach model was formulated in the paper^3Sl f or analysis of the role played by 
layered structures in superconducting materials IM2H342| -p^e mo del considers a stack of parallel 
two dimensional superconducting layers separated by an insulated material (or vacuum), with 
a nonlinear interaction between the layers. It is also assumed that an external magnetic field 
is applied to the system. In some sense, the Lawrence-Doniach model can be considered as 
an anisotropic version of the Ginzburg- Landau model) 31 ™ 338 ! More specifically, an anisotropic 
Ginzburg-Landau model can be considered as a continuous limit approximation to the Lawrence- 
Doniach model. However, when the coherence length in the direction perpendicular to the layers is 
less than the distance between the layers, these models are difficult to compare. In the framework 
of the approach used by Lawrence and Doniach the superconducting properties of the layered 
structure were considered under the assumption that in the superconducting state the free energy 
per cell relative to its value in the zero external field can be written in the following form 



m = E 



a l {T)\* i (r)\'+P\^ l {r)\ A + 



1 



2m ab 



iftY + — ]*i(f) 



<ij> 



1.2) 



Here, ^(r) is the order parameter of the Ginzburg-Landau order of the layer number i, (*$>i(x,y) 
is a function of two variables), the operator V acts in the x-y plane; A is the corresponding 
vector's potential, a and f3 are the usual Ginzburg-Landau parameters, rjij describes a positive 
Josephson interaction between the layers; and < ij > denotes summation over neighboring layers. 
It is assumed that the layers correspond to planes ab, and the c axis is perpendicular to these 
planes. Accordingly, the z axis is aligned with c, and the coordinates x-y belong to the plane ab. 
The quantities rjij are usually written as follows 



h 2 



2m r s 2 



(8.3) 



Here, s is the distance between the layers. As one can see, for a rigorous treatment of the 
problem one has to take into account the anisotropy of the effective mass at the planes ab and 
between them, m a i, and m c , respectively. Frequently, the distinction between these two types of 
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anisotropy is ignored, and a quasi-isotropic case is considered. If we write down \&j in the form 
= \^fi \ exp(i(fi) and assume that all \^S?i\ are equal, then rjij is given by 

h 2 

= O 2 l^l^ 1 ~ COS (^ ~ Vi-l)]" ( 8 ' 4 ) 

The coefficient a«(T) for the layer number i is given by 

a t (T) = ~P , (8.5) 

i 

where Tf denotes the critical temperature for the layer number i. Next, one can consider the 
situation where VPj(r) = VPi(r) and A = 0. In the vicinity of T c the contribution from is 
small. Taking into account all these simplifications one can write down the free energy's density 
in the following form 

n 

f = Y,a l (T)\* i \ 2 + ^ //;,•*; V- (8.6) 

i <ij> 

This is the quasi-isotropic approximation with single mass parameter a. The Ginzburg-Landau 
equations follow from the free-energy extremum conditions with respect to variations of 



6% 



(ati + rji-i i + rji - {rn-i + rn = 0. (8.7) 



The corresponding secular equation is given by 



(aj(T) + rji-i i + rji i+i)<% - rjij5i j± i 



0. 



It is assumed in the framework of the Lawrence-Doniach modeP^ that the transition temperature 
corresponds to the largest root of the secular equation. In other words, one has to investigate 
solutions of the equation 



a i u i 



0, (8.9) 



or, in other form 



det(T/ — M) = 0, where M i5 = (if - ^=^Tf - ^^T t ° ) <5« + ^lfS ij±1 . (8.10) 

Thus, the problem is reduced to finding the maximal eigenvalue of the matrix M. If we take into 
account the external field, then the complete form of the Lawrence-Doniach equatiorP^J is given 
by 

ti 2 ( 2e A 2 

a% + P\*i\ 2 *i-- (v + i—A) % 

2m ab V nc I 



ti 



2 



2m r s 2 



A large number of papers are devoted to investigations of the Lawrence-Doniach model and 
to development of various methods for its solution .1MQH345J j n mari y respects this model corre- 
sponds to layered structures of high-temperature superconductors)*^ and in particular to mer- 
curocuprates I342H342I ^ relativistic version of the Lawrence-Doniach model was studied in the pa- 
p er |3i3| wriere violation of the local U{\) gauge's symmetry was considered by analogy with Higgs 
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mechanism PS3 A spontaneous breaking of the global U(l) invariance is taking place through the 
superconducting condensate. The paper^21 a l so studied in detail the conse que nces of spontaneous 
symmetry breaking in connection with the Anderson-Higgs phenomenonP^I As was mentioned 
already, the concept of spontaneous symmetry breaking corresponds to situations with symmetric 
action, but asymmetric realization (the vacuum condensate) in the low-energy regime. As a result 
the realization has a lower symmetry than the causing actionPGSEEI] 

In essence, the Higgs mechanisrrPE! follows from the Anderson idecP^J Q n the connection between 
the gauge's invariance breaking and appearance of the zero-mass collective mode in superconduc- 
tors. Difference-differential equations for the order parameter, as well as for the vector potential at 
the plane and between the planes were also derived in the paperF^ These equations correspond 
to the Klein-Gordon, Proca and sine-Gordon equations. The paper also contains a comparison 
of the superconducting phase shift {ipi — '■Pi-i) between the layers in the London limit with the 
standard sine-Gordon equation. A possible application of this approach to description of the 
high-temperature superconductivity in layered cuprates with a single plane in the elementary cell 
and with a weak Josephson interaction between the layers was also considered. 
Thus, a systematic scheme for a phenomenological description of the macroscopic behavior of 
layered superconductors can be constructed by applying the covariance and gauge invariance 
principles to a four-dimensional generalization of the Lawrence-Doniach model. The Higgs mech- 
anism®^ plays the role of a guiding idea, which allows one to place this approach on a deep and 
nontrivial foundation. The surprising formal simplicity of the Lawrence-Doniach model once again 
stresses the R. Peierls idedP^ on the efficiency of physical model creating. 



9 Nonequilibrium Statistical Operators and Quasiaverages in the 
Theory of Irreversible Processes 

It has been mentioned above that Bogoliubov's quasiaverages concepP^"^' plays an important 
role in equilibrium statistical mechanics. According to that concept, infinitely small perturba- 
tions can trigger macroscopic responses in the system if they break some symmetry and remove 
the related degeneracy (or quasidegeneracy) of the equilibrium state. As a result, they can pro- 
duce macroscopic effects even when the perturbation magnitude is ten d to ze ro, provided that 
happens after passing to the thermodynamic limit. D.N. Zubarev showed 3 -^"*^ 8 -' that the concepts 
of symmetry breaking perturbations and quasiaverages play an important role in the theory of 
irreversible processes as wellP^ The method of the construction of a nonequilibrium statistical 
operator^ becomes especially deep and transparent when it is applied in the framework of the 
quasiaverage concept. The main idea of the papers^^"^ 8 -' was to consider infinitesimally small 
sources breaking the time-reversal symmetry of the Liouville equation 

5eM + l w ,o) lfl i = o <,.,, 

which become vanishingly small after a thermodynamic limiting transition. 

The main idea of the method of a nonequilibrium statistical operator (NESO) 38 can be summarized 
as follows. In the scale of sufficiently large times the nonequilibrium state of the system can be 
described by some set of parameters F m (t), and one can find such a particular solution of the 
Liouville equation (19. ID which depends on time only through F m (t). The first argument of the 
operator p{t, 0) refers to an implicit time dependence. It is assumed that the nonequilibrium 
statistical ensemble can be characterized by a small set of relevant operators P m {t) (quasi- 
integrals of motion). The corresponding NESO is a functional of P m {t): 

p(t) = p{...P m (t)...}. (9.2) 
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One can show, seej^that if the statistical operator p(t, 0) satisfies the Liouville equation, then it 
is given by 

p = exp(A- / dti J^G m (ti)P m (ti)); A = 1 - A, (9.3) 

J-co m 

where 

G m (t 1 )=ee^F m (t + t 1 ), (9.4) 

A = e dt ie £tl X(t + ti) = A(t) - / dtie Etl A(i + ti). (9.5) 

J —CO J —OO 

Alternatively, it can be written as follows 

/o iHt iHt 

dt 1 e £tl e ( —^)lnpg(t + t 1 )e ( - — -^)) = (9.6) 

fO 

exp(-S(t, 0)) = exp(-£ / dtie etl S'(i + = 

J — OO 



-OO 

r0 



exp(-5(t,0)+ / dtie^^t + ti,*!)), 

J —CO 



where 



Pq = exp I n-J2 F m(t)P m J =exp(-5(i,0)), O = InTr exp I - ^ F m (t)P m j 
\ m / \ m / 



(9.7) 



5(i,0) = ^M+ 1 [S ( t ,0), H]; (9.8) 
at zn 

S(t, tl ) = exp( i -^) S(t, 0) exp ( ~ ' '' ' ' 1 



h J w ' r V 

Here, /9 q is the quasi-equilibrium statistical operator, which corresponds to the extremum value 
of the information's entropy 

S = -Tr(plnp), (9.9) 
under the additional conditions that Tr(pP m ) = (P m )q', Trp = 1. In this case 

8$ 

-(P m ) q ; (...)« = T±(pa...), (9.10) 



SF m 

$(p) = -Tr(plnp) - ^F m Tr(pP m ) + ATrp, (9.11) 

m 

(Pm)' = {Pm)\. (9-12) 

The quantum Liouville equation (|9.ip (as well as the classical one) is invariant with respect to the 
time reversal. One can show^HGim that p(t, 0) satisfies the Liouville equation with an additional 
infinitesimally small (proportional to s) source-term in the right hand side, and we send e to zero 
after the thermodynamic limit. Indeed, let us consider the equation 

^l + l[p £ ,H] = -e(p £ - Pq ) (9.13) 

or, equivalently, 

+ ^[lnp £ ,H] = -e{lnp e - hxp q \ (9.14) 
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where e — > after passage to the thermodynamic limit. Equation (|9.13p is an analogue of the 
corresponding equation in the quantum scattering theory GSSlElQ] The introduction of infinitely 
small sources in the Liouville equation corresponds to imposing the following boundary conditions 

( iHt\\ , , . . ,. / —iHt\\ 
exp ( J (p(i + ii)-p g (t + ti))exp ( — ^ J 0. (9.15) 

Here,ii — > — oo after the thermodynamic limit. It was shown in the papers^^^IKEI that the 
operator/},, is given by 

p £ (t,t)=e dtie^ 1 -^ p q (tx,t x ) = £ dhe £tl p q {t + t u t + h). (9.16) 

J —oo J —oo 

Here, the first argument in p(t, t) refers to the implicit time dependence via the parameters F m (t), 
while the second argument refers to the time dependence via the Heisenberg representation. The 
desired statistical operator is given by 



f0 

p £ = p e (t,0) =p q (t,0) =e / dtxe^Pqit + hM)- 

J —oo 



(9.17) 



(9.18) 



Hence, the nonequilibrium statistical operator is given by 

p = Q~ 1 exp \-J2 B ™ J =Q" 1 exp (-J^e / dt ie £tl (F m (t + h)P m (ti)) J = 

Q~ 1 exp(-^2F m (t)P m + '^2 j dhe £tl [F m (t + *i)P m (ii) + F m (t + *i)P m (*i)] J • 

One can rewrite Eq. (|9.14p in the following form 

— (e et In p(t,t)) = ee et In p q (t,t), (9.19) 

where 

\np(t,t) = U\t,0)\np(t,0)U(t,0); U(t, 0) = exp V (9.20) 
Integrating Eq. fj9. 191) over the interval (— oo,0) we obtain 

r° 

]np(t,t)=e dtie £tl \np q {t + ti,t + ti). (9.21) 

J —oo 

It is assumed that lim e ^. + ^ n p(t> t) = 0- Therefore, 

p(t,0) = exp (- £ J dhe 611 lnp q (t + t 1 ,t 1 ) \ = exp (lnp q (t, 0)) = exp {-S(t, 0)). (9.22) 

The average value of any dynamic variable ^4 is now given by 

(A) = lim Tr(p(t,0U). (9.23) 

£->0+ 

We see that the above average is in fact nothing else but a quasiaverage. The normalization of 
the quasi-equilibrium distribution p q is preserved if 

Te(p(t, 0)P m ) = (P m ) = (P m ) q ; Tip = 1. (9.24) 

Thus, one can assert that the origin of the irreversibility effect is closely related to the viola- 
tion of the time-reversal symmetryp^U as well as to the notion of quasiaverages from statistical 
mechanics EHElZlEa 
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9.1 Generalized Kinetic Equations 

The NESO method 38 found wide applications in various problems of statistical mechanics. An 
important contribution to the development of the kinetic equations' theory in the framework of 
NESO method was made by L.A. Pokrovskii I352H354I Generalized kinetic and transport equations 
describing the time evolution of the variables (P m ) and F m (t) are obtained by averaging the 
equation of motion for P m over the derived NESO 

The generalized transport equations are given by 

The corresponding entropy production can be written down in the following form 

S(t) = (S(t,0)) =-J2(Pm)F m (t) = ~E af mL (f M t)F m (t). (9.27) 

n,m or m^)Or n yi) 

The two equations in (|9.26p are mutually conjugate, and together with Eq. (|9.27p they form a 
complete system of equations for calculation of the quantities (P m ) an d F m . 
Now, following the paper j 3 -^ we are going to write down kinetic equations for a system with weak 
interaction. The corresponding Hamiltonian is given by 

H = Hq + V. (9.28) 

Here, Hq is the Hamiltonian of noninteracting particles (or quasiparticles) and V is the interaction 
operator. As the set of relevant operators we choose the operators P m = Pk of the form a k a k or 
a\.a k+q . Here, at and a k are the usual creation and the annihilation operators (either Fermi or 
Bose). We begin with the following equations of motion: 

P k = hp kl H]. (9.29) 

in 

It is usually assumed that 

[P k ,H ] = J2 c ki p U (9.30) 
i 

where c k i are some coefficients (c- numbers). 
According to Eq. (|9. 18|) we have 

p = Q~ 1 exp l-Y^F k (t)P k + Y, I dhe^ihit + t^P^ + Fkit + t^Pkih)}) . (9.31) 
V k k J ~°° J 



Keeping in mind that (P k ) = (P k ) q , we can write down the generalized kinetic equations^' for 
(P k ) as follows 

The right hand side of Eq. (|9.32|) contains the generalized collision integral, which, using an 
expansion in powers of V, can be written as follows 

^ = Ll + Ll + Lf + Lf, (9.33) 



52 



where 

i 

Ll = ±([P k ,V]) q , (9.35) 

Lf = ^ J^dhe^dVih),^)])^ (9.36) 

Analogously one can find the higher-order terms V 3 , V 4 and so on. 



9.2 Generalized Kinetic Equations for a System in a Thermal Bath 

The papers^SHlSZ] ^ see gjgjlt) generalize the equations in (|9.32p for the case of a system interacting 
with a thermal bath. The concept of a thermal bath or heat reservoir is fairly complicated and 
has certain specific featuresP^ According to the standard definition, a thermal bath is a system 
with, effectively, an infinite number of degrees of freedom. A thermal bath is a heat reservoir 
maintaining the investigated system under a particular temperature. Following Bogoliubovp^ we 
will assume that a thermal bath is a source of stochasticity for a small subsystem (which in an 
extreme situation can be just a single particle). Such a small subsystem can be, for example, an 
atomic or a molecular system interacting with an electromagnetic field, or a system of nuclear or 
electron spins interacting with the crystal lattice, etc. We will describe the entire system by the 
Hamiltonian 

H = H 1 + H 2 + V, (9.38) 

where 

Hx = J2 E«a\fla, y = Yl 'bapatap, $ Q(3 = $]j a . (9.39) 

a a,/3 

Here, Hi is the Hamiltonian of the small subsystem; the creation and annihilation 

operators of quasiparticles with the energies E a in the small subsystem; V is the operator describ- 
ing the interaction between the small subsystem and the thermal bath; and H 2 is the thermal 
bath's Hamiltonian, which W6 do not write down explicitly. The quantities *&q>b 

are operators 

acting on the thermal bath's degrees of freedom. We assume that the state of the system can 
be characterized by a set of operators (P a p) = {a a ap), and the state of the thermal bath by the 
operator {Hq). Here, (. . .) denotes the averaging with respect to the NESO, which is defined as 
follows: 

p g (t) = exp(-S(t,0)), S(t,0) = n(t) + Y,PapF a p(t) + !3H 2 , (9.40) 

a/3 

O = InTr exp(- P a pF a p{t) - j3H 2 ). 

af3 

Here, F a p{t) are the thermodynamic parameters conjugate to P a p, f3 is the inverse temperature of 
the thermal bath. All operators are considered in the Heisenberg representation. We write down 
the nonequilibrium statistical operator as follows 



p(t)=exp(-S(t,0)), (9.41) 
5p) = e ^ dt ie £t i ( fl(t + t x ) + P a pF a p{t) + (3H 2 
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The parameters F a p(t) are determined by the condition (P a p) = {P a /3)q- To derive the kinetic 

equations we will use an expansion over the small parameter in the interaction V. It is also 
assumed that the equation {& a /3)q = 0. holds. It is convenient to rewrite p q as follows 

p q = P1P2 = Qq 1 exp(-L (t)), (9.42) 

where 

Pi = Qi l exp I - P a pFaf)(t) ; Qi = Tr exp I - ^ P a pF a p(t) , (9.43) 

\ a/3 / \ a/3 / 

P2 = Q2 1 e~ PH2 ] Q 2 = Tiexp(-/3H 2 ), (9.44) 

Q q = QiQ 2 ] L Q = J2 p apFap(t)+pH 2 . (9.45) 

a/3 

We begin from the following relationship: 

= ^([P a p,H]) = -(^ - E a )(P aP ) + ^{[P^V]). (9.46) 

We terminate the expansions at the second order terms in V . The kinetic equations for the 
quantities (Pa/3) of a system in a thermal bath are given by 



d(Pgp) 

dt 



i i r° 

-{E p - E a )(P aP ) - J dtu 



(\[P a p,V\,V{h)\) q . (9.47) 



These equations generalize the results of the paper^^l f or a system in a thermal bath. One can 
show that the choice of the concrete model's form for the Hamiltonian (|9.38[) is not essential. For 
arbitrary H\ and V, and for some set of variables (P&) satisfying the condition [Hi, P^] = Y2i c klPu 
one can construct a quasi-equilibrium statistical operator p q in the following form 



Pg 



Qq 1 exp (- PkFkff) - (3H 2 ^ . (9.48) 



Here, Fp.(t) are the parameters conjugate to (Pfc). The kinetic equations for (P^) are given by 
d(Pk) ' 1 r0 



dt 



lJ2<*i(Pl)-]pJ_ dhe^m^Vih)]),. (9.49) 



9.3 A Schrodinger-Type Equation for a Dynamic System in a Thermal Bath 

Following the papersPEHJESl we consider now the behavior of a small dynamical subsystem with 
a Hamiltonian H\, which interacts with a thermal bath described by the Hamiltonian H 2 . As 
the operators characterizing the state of the small subsystem we choose the operators a a ,a a , and 
n a = a a a a . In this case the quasi-equilibrium statistical operator p q is given by 



p q = exp (n - Y,(f*(t)a a + + F a (t)n a ) - (3H 2 ^ = exp(-5(t, 0)), 

Q = InTrexp ^- ^(/ a (i)a Q + + F a (t)n a ) - fiH^j 



(9.50) 
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Here, f a , f a and F a , play the role of Lagrange multipliers. They are the parameters conjugate to 
(a a ) q , (a a ) q and (n a ) q : 

i \ ^ i \ ^ aa ^ tp i ,„ 

[O'a/q = ~~ x j, 7TT) \ n a)q = — e j-i 77\"> ~Fj \~ = Ja(t), T7 \~ = ^a\t)- (9.51 J 

ofa{t) oF a (t) d(a a ) q o{n a )q 

The quantities a a , a a in the statistical operator can be interpreted as sources of quantum noise 
(those terms break the spin conservation law in the case of fermions; see the papers^^) . Let us 
write down the quasi-equilibrium statistical operator as follows 

Pq = P1P2, (9.52) 

where 



(9.53) 



Pi = exp - £}(/ a (i)a a + fl(t)ai + F a (t)n a )j 

fii = InTrexp ^- ^(/ a (t)a a + fl{t)a) a + F a (t)n a )\ 

p 2 = exp (0 2 - pH 2 ) , fi 2 = InTrexp (-/3H 2 ) (9.54) 

As a result we obtain the expression (I9.4ip for the NESO p. We assume that the following 
conditions are satisfied: 

(a a ) g = (a a ), (a a ) g = (a a ), (n a ) q = {n a ). (9.55) 
We start from the equations of motion 

ih^- = ([a a ,H 1 ]) + ([a a ,V]), (9.56) 

In the second order in V we obtain 

_d{a, 



in- 



dt 

_d(n 



E a (a a ) + ^-[ dhe^iiia^VlVih^q, (9.58) 
^ y° ^ ie rfl ([[n Q ,F],y(ti)]) 9 . (9.59) 



in 



dt 

Here, V(t\) denotes the operator V in the interaction representation. The expansion yields 

d(a a ) 
iri—jj— = E a (a a ) + 

-^J dtie etl I ^{Qap^ixAti^qiapala^q - (^y, v {t{)^ a p) q (a^a v ap) q J , (9.60) 
where ^(ti) = exp(|(^ - E v )ti), or, equivalently 



1 f 

j E a (a a ) + -J2 dtie etl (^(t> fl p(t 1 ))q{a^+ (9.61) 

1 f° 

/ ^ie e * 1 ([# ai ,,0^(ti)])q(aj e a i ,a^)g. 
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Therefore, we obtain 

d(a a ) 



ih- 



1 f° 

E a {a a ) + - i J2 dhe^i^a^ih^qiap). (9.62) 

Using the spectral representations for the correlation functions we can write down 

ifggt = Ea (a a ) +J2K af) (ap), (9.63) 
at 

where K a p are defined as follows: 



^X)/ *ie dl <^V(*i)>?= ( 9 - 64 ) 

9^ ^ 



It 



Thus, we have obtained a Schrddinger-type equation for the mean amplitudes {a a ). In a certain 
sense this equation is an analogue (or a generalization) of the Schrodinger equation for the case 
of a particle moving in a medium. Let us consider this analogy in a more detail. First, we write 
down the analogue of the wave function in the following form 

W) = Y,X a (r f )(a a ). (9.65) 

a 

Here, {Xa(^)} is a complete orthonormal set of single-particle eigenfunctions of the operator 
( — ^m^ 2 + v (^)j > where v(r) is the potential energy, 

h 2 \ 

-— V 2 + v{f)\ Xa (f) = E aXa {f). (9.66) 

Thus, the quantity ip(r) plays the role of a wave function describing a particle moving in a medium. 
Equation (|9,63p can be rewritten in the following form 

.^a^r) = f_^_ v 2 + y , ff) \ + f K(?,pU(?)dP. (9.67) 
at \ 2m J J 

The kernel K(r, r') of the integral equation (|9.67p is given by 

K(?y) =J2K a /3X<*(7)4(r') = 4 E T *i^<*«/^(*i)>gX«Wx^0- (9.68) 

We see that Eq. (|9.67p can indeed be classified as a Schrodinger-type equation for a dynamical 
system in a thermal bath. It is interesting to note that very similar equations of the Schrodinger- 
type with a nonlocal interaction were used in the collision theorySSSl for description of particle 
scattering on a cluster of many scattering centers. 

In order to make clear some special features of Eq. (|9.67|) let us consider the translation operator 
exp(iqp/H), q = r' — r, where p = —iKV r , Then, Eq. (|9,67p can be rewritten in the following form 

.^r) = ^_|1 V 2 + v( ^\ ^ + £ Dfrfi^ (9.69) 
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where 

d 3 qK(r,r + q)e—. (9.70) 

It is reasonable to assume that the wave function ip(r) does not change very rapidly over distances 
comparable to the characteristic correlation length of the kernel K(r,r'). Then, using the series 
expansion for exp(iqp/h) in Eq. ([9.69P we obtain in the zeroth order 

if d^r)_ = f_^_ v 2 +v , f * + R e u(r)]^(r) + ilmU(rU(r), (9.71) 
ot \ 2m J 

where 

U(r) = Re U(r) + ilm U(f) = j SqK(r,r + q). (9.72) 

Equation (|9.7ip has the exact functional form of a Schrodinger equation with a complex potential 
well known in the collision theory P^S! Note that the introduction of the quantity ip{r) does not 
mean that the state of the small dynamical subsystem becomes pure. The state remains mixed 
because it is described by a statistical operator. The dynamics of the system is described by a 
system of coupled evolution equations for the quantities f a , f a and F a . Note, that there were many 
attempts to derive a Schrodinger- type equation for a particle in a medium 1261H363J KorringcP^ tried 
to obtain such an equation in the form of an evolution equation with a non-hermitian Hamiltonian. 
However, his equation (cf. Eq. (29) fronP^) 

^ = (*w+W + if...)r Wl ( 9 .73, 

where W'{t) is a statistical density matrix describing the original system, is rather a modified Bloch 
equation. An attempt to derive a Schroedinger-type equation for a Brownian particle interacting 
with a thermal environment was made in the paperP^J Tli e evolution equation obtained there is 
given by 

where 

W(t) = -(|£) J r H^Wdr. (9.75) 

Here, / is the friction coefficient, Vr is a random potential, and Vn(r, t) = —rFn(t), where Fji(t) is 
a random vector- function of time. Excluding the function W(t) with the help of the transformation 

ij)(r, t) = exp[i9(t)](t){f, t), (9.76) 

where 

0(f) = _^- 1 exp(-^) / eM^-)W(s)ds, 
m J t m 

one obtains the equation for (p(r, t) in the following form 

(r,f)m[- 



2im ' (f)*(f,t) 

It is clear that the dynamic behavior of a particle in a dissipative environment is most accurately 
described by a Schrodinger- type equation with damping, see Eq. (|9.67|) above. This is actually 
the reason for applications of this equation in numerous problems of physics, physical chemistry, 
biophysics, and other areas JHHIH3Z2I a more detailed discussion of various aspects of dissipative 
behavior and of stochastic process in complex systems is given in the reviews J3I222H22SI 
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10 Conclusion 



In the present paper we have shown the determining role played by correlation effects in system- 
atic microscopic descriptions of magnetic, electrical, and superconducting properties of complex 
substances. We have stressed that the approximation of tight-binding electrons and the method 
of model Hamiltonians are very effective tools for description of these substances. In many cases 
the methods of quantum statistical mechanics, many of which were formulated and developed by 
N.N. Bogoliubov, allow one to develop efficient approaches for solution of complicated problems 
from microscopic theory of correlation effects, especially in the case of strong electron correlations. 
The method of two-time temperature Green's functions allows one to efficiently investigate the 
quasiparticle dynamics generated by the main model Hamiltonians from the quantum solid state 
theory and the quantum theory of magnetism. The method of quasiaverages allows one to take 
a deeper look at the problems of spontaneous symmetry breaking, as well as at the problems 
of symmetry and dissymmetry in the physics of condensed matter. Further development of the 
theory describing many-particle effects and investigations of more realistic models will allow one 
to gain more precise ideas on the effective interactions in the systems, which determine various 
phenomena, the main features of electron states, and therefore, the physical properties of real 
substances. The methods developed by N.N. Bogoliubov are and will remain the important core 
of a theoretician's toolbox, and of the ideological basis behind this development. 
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